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ABSTRACT 

We analyze 221 eclipsing binaries (EBs) in the Large Magellanic Cloud with B-type main-sequence 
(MS) primaries (Mi ss 4-14 Mg) and orbital periods P = 20-50 days that were photometrically 
monitored by the Optical Gravitational Lensing Experiment. We utilize our three-stage automated 
pipeline to (1) classify all 221 EBs, (2) fit physical models to the light curves of 130 detached 
well-defined EBs from which unique parameters can be determined, and (3) recover the intrinsic 
binary statistics by correcting for selection effects. We uncover two statistically significant trends with 
age. First, younger EBs tend to reside in dustier environments with larger photometric extinctions, 
an empirical relation that can be implemented when modeling stellar populations. Second, younger 
EBs generally have large eccentricities. This demonstrates that massive binaries at moderate orbital 
periods are born with a Maxwellian “thermal” orbital velocity distribution, which indicates they 
formed via dynamical interactions. In addition, the age-eccentricity anticorrelation provides a direct 
constraint for tidal evolution in highly eccentric binaries containing hot MS stars with radiative 
envelopes. The intrinsic fraction of B-type MS stars with stellar companions q = M 2 /M 1 > 0.2 
and orbital periods P = 20-50 days is (7 ±2)%. We find early-type binaries at P = 20-50 days are 
weighted significantly toward small mass ratios q « 0.2-0.3, which is different than the results from 
previous observations of closer binaries with P < 20 days. This indicates that early-type binaries 
at slightly wider orbital separations have experienced substantially less competitive accretion and 
coevolution during their formation in the circumbinary disk. 

Subject headings: binaries: eclipsing, close; stars: massive, formation, evolution, statistics 


1. INTRODUCTION 

It has long been understood that the main-sequence 
(MS) binary star fraction increases with primary mass 
(lAbtjll983l: iDuauennov fc M a vail I 199 ll: iFischer fc Marcvl 
Il992t iRaghavan et al l 120101 iDuchene fc Kraus! 120131 
etc.). Indeed, most massive stars with M\ > 
10 M 0 will interact with a stellar companion before 
t hey explode as core-collapse supernovae (IS ana et al.1 
I2012D . Throughout the decades, there have been 
significant advances in the detection of clo se and 


wide companions to massive stars (|Wolff 

1978; 

Garmanv et aL| 1980| Levato et al. 19871: Abt et al. 

1990: 

Shatskv & Tokovinin 2002| Kouwenhoven et al. 

2007; 

Sana et al. 201(1; Rizzuto et al. 20131; Kobulnickv 

et al. 

2014[). However, the intrinsic properties of binary 


companions to early-type primaries, e.g. their 
eccentricity and mass-ratio distributions, remain elusive 
at intermediate orbital periods. The major goal of 
this work is to help fill this particular portion of the 
parameter space. 

Eclipsing binaries (EBs) offer a key to the accurate 
measurement of the binary properties of early-type stars. 
Large photometric surveys, such as the third phase of the 
Optical Gravitational Lensing Experiment ( OGLE-I II), 
have discovered tens of thousan ds of EBs (|Graczvk et al.1 
120111 IPietrukowicz et all 120131) . These populations of 
EBs are orders of magnitude larger than previous binary 
samples. Despite the geometrical selection effects, we can 
still achieve large sample statistics to reliably infer the 
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intrinsic binary fraction and properties at intermediate 
orbital periods. We emphasize that EBs can probe a 
unique portion of the binary parameter space unavailable 
to ot her observational techniqu es. 

In IMoe fe Di Stefanol (|2013l hereafter Paper I), we 
incorporated OGLE catalogs of EBs in the Large and 
Small Magellanic Clouds (LMC and SMC, respectively) 
as well as Hipparcos observations of EBs in the 
Milky Way. We compared the close binary properties 
{P < 20 days) of early-B MS primaries in the three 
different galaxies. The Milky Way and SMC EB samples 
are too small to warrant an analysis of period-dependent 
binary propert ies. The OGLE-III LMC EB catalog 
(jGraczvk et all 120111) . on the other hand, contains 

5 - 40 times more systems, is relatively complete toward 
shallow eclipse depths, and includes the full I-band and 
V-band ligh t curves. _ 

In IMoe fc Di Stefanol (|20l4 hereafter Paper II) , we 
developed a three-stage automated pipeline to analyze 
EBs with short orbital periods in the OGLE-III LMC 
database. This pipeline (1) classifies EBs according 
to their light curve characteristics, (2) measures the 
intrinsic physical properties of detached EBs, e.g. ages 
and component masses, based on the observed radii 
and temperatures, and (3) recovers the intrinsic binary 
statistics by correcting for selection effects. 

In the present study, we utilize EBs in the OGLE-III 
LMC database to measure the binary fraction, mass-ratio 
distribution, and eccentricity distribution of B-type MS 
stars with intermediate orbital periods P = 20-50 days. 
We organize the rest of this paper as follows. In §2, 
we define our selection criteria for identifying EBs 
with B-type MS primaries, intermediate orbital periods, 
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and well-defined eclipse parameters. We next describe 
an automated procedure we developed to fit detailed 
physical models to the observed EB light curves, and 
we present our results for the physical properties of 
the individual EBs (§3). In §4, we explain the 
observed trends in the measured EB parameters, paying 
special attention to the empirical age-extinction and 
age-eccentricity anticorrelations. We then perform 
Monte Carlo simulations to quantify selection effects, and 
present our results for the corrected binary statistics (§5). 
We summarize our main results and conclusions in §6. 

2. EB SELECTION AND CLASSIFICATION (STAGE I) 

In Paper II, we developed a three-stage automated 
pipeline to fully analyze short-period EBs in the 
OGLE-III LMC database. In the present study, we 
adapt our routine to identify intermediate-period EBs 
with well-defined light curves (Stage I - this section), 
measure their physical properties (Stage II - §3), and 
correct for selection effects (Stage III - §5). EBs 
with intermediate orbital periods exhibit two major 
differences that must be considered. First, the eclipse 
widths 0i and 0 2 , which are expressed as a fraction of 
the orbital period P, become narrower with increasing 
orbital separation. Given the average number (A/)) 
~ 470 of I-band measurements in the OGLE-III LMC 
survey l|Graczvk et all 1201111 , the light curves are not 
sufficiently sampled if either of the eclipse widths 
0 < (A//) -1 « 0.0021 are too narrow. EBs with small 
MS components and long orbital periods P > 50 days 
have narrow eclipses 0 < 0.002, and are therefore 
not Nyquist sampled. This subsampling leads to 
detection incompleteness, issues with aliasing, and the 
inability to fully characterize their intrinsic physical 
properties. Hence, it is the finite cadence of the 
OGLE-III observations, not geometrical selection effects, 
that limi ts our present stu dy of EBs to P = 20-50 days 
(see also lSbderhielm fe Dischleill2005lh 

Second, the majority of early-type EBs at P > 20 days 
are in eccentric orbits. We must therefore adapt our 
physical models to simultaneously fit the eccentricity e 
and argument of periastron u (§3). In addition, 
it is possible for an eccentric binary to have a 
certain combination of eccentricity, periastron angle, 
and inclination that is sufficiently offset from edge-on 
(e.g., * < 86°) so that there is only one eclipse per orbit. 
Indeed, there are many EBs with single eclipses in the 
OGLE-III LMC database (see below). Unfortunately, 
we cannot measure the physical properties of these 
systems. We therefore remove single-eclipse EBs from 
our well-defined sample, and we account for their removal 
when we correct for selection effects (see §5). In the 
following, we review our methods from Paper II, where 
we pay special attention to the nuances of EBs with 
intermediate orbital periods. 

In this study, we select the A/r ~ 96,000 syste ms in 
the OGLE-III LMC catalog (jUdalski et all l2008f) with 
mean magnitudes 16.0 < (I) < 17.6 and observed colors 
—0.25 < (V — I) < 0.20. Given the distan ce modulus 
fj = 18.5 to the LMC (jPietrzvnski et alJ2013f l and typical 
dust reddenings E(V — I) ~ 0.1-0.3 mag tow ard hot 
young stars in the LMC (jZaritskv et al.l 1200411 . these 
stars have luminosities and surface temperatures that 
correspond to B-type MS primaries. From this sample, 


we analyze the 221 systems that were id entifie d as EB s 
with orbital periods P = 20-50 days (|Graczvk et al.l 
|20H . In Table 1, we list the OGLE-III LMC EB 
identification numbers, observed colors (V — /), and 
numbers of I-band measurements A// for each of these 
221 EBs. 

As in Paper II, we measure the intrinsic rms scatter in 
the /-band light curve outside of eclipses for each EB. We 
then calculate the correction factor f a j > 1.0, i.e. the 
ratio between the actual rms scatter and photometric 
uncertainties reported in the catalog. For each /-band 
measurement in an EB light curve, we multiply the listed 
photometric uncertainties by the correction factor f a j to 
determine the corrected uncertainties. 

We classify EBs based on an analytic light curve model 
of two Gaussians with eight total free parameters. The 
orbital phase 0 < <fi < 1 is determined by the time 
of observation and two model parameters: the orbital 
period P (in days) and epoch of primary eclipse minimum 
t a (Julian date — 2450000). The six remaining analytic 
model parameters are the average I-band magnitude 
outside of eclipses (/), primary and secondary eclipse 
depths All and A/ 2 , primary and secondary eclipse 
widths 0i and 0 2 , and the phase of secondary eclipse 
<f> 2 . The analytic model of Gaussians is: 


IgW 


1 2 

(I) + Ah ex p( 2 q 2 ) + exp( 


A/ 2 exp ^ 


-(</> - <f> 2 ) s 
202 


~( 0 ~ 1) 2 

202 


(1) 


We fit this analytic model to each EB I-band 
light curve. Specifically, we utilize an automated 
Levenberg-Marquardt algorithm ( MPFIT, iMarkwardtl 
120091) to minimize the Xg statistic. The MPFIT routine 
provides robust best-fit solutions and measurement 
uncertainties for the eight analytic model parameters. 
Some of the photometric measurements are clear outliers, 
so we clip up to A/" c < 2 data points per light curve 
that exceed 4<r from the model. This results in 
v = Mi — M c — 8 degrees of freedom. For each EB, 
we report in Table 1 the eight fitted analytic model 
parameters and the fit statistics. Excluding the few EBs 
that exhibit variability or are evolved Roche-lobe filling 
systems (see below), the goodness-of-fit statistics Xg/ 11 
= 0.87 -1.16 indicate the analytic models can adequately 
describe the EB light curves. 

We can measure the physical properties of EBs based 
solely on the observed photometric light curves (see 
§3) only if: (1) the binary components are detached 
from their Roche lobes, (2) the light curves have two 
well-defined eclipses, and (3) there is no superimposed 
variability. To be considered well-defined, we require 
that the ler uncertainties in the measured eclipse depths 
Ah and A/ 2 and eclipse widths 0i and 0 2 are <20% 
their respective values. These criteria are not satisfied 
for 91 of the 221 EBs due to a variety of reasons, which 
we discuss below: 

(A) No Secondary Eclipse. For 16 of our EBs, there is 
no evidence for a secondary eclipse. These EBs may have 
secondary eclipses that are too shallow and below the 
sensitivity of the OGLE-III LMC survey, or have eclipse 
widths that are too narrow and therefore not detected 
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Fig. 1.— Six examples of the 91 EBs that have properties that are uncertain, variable, peculiar, etc., which leaves 130 EBs in our 
well-defined sample. Panel A: one of 16 EBs that does not have a visible secondary eclipse. Panel B: one of 32 EBs that show both eclipses 
but where one of them is too narrow and/or too shallow to be accurately measured. Panel C: one of the three EBs with wide eclipses 
that demonstrate one or both components fill their Roche lobes. Panel D: one of the 23 EBs with an ambiguous orbital period. Using the 
catalog orbital periods (black), these systems have nearly identical eclipses separated by almost precisely 50% in orbital phase. The more 
plausible scenario is that these EBs have half the listed orbital periods (red) and therefore exhibit one eclipse per orbit such as the example 
shown in panel A. Panel E: one of the 15 EBs that exhibit variability. Three of these systems are intrinsic variables. The other 12, such 
as the displayed example, show changes in the eclipse properties most likely caused by orbital motion with a tertiary companion. Panel F: 
one of the two EBs with peculiar light curve properties. 


given the cadence of the observations. Most likely, the 
EBs have a certain combination of e, ui, and i as discussed 
above so that there is only one eclipse per orbit. We list 
these 16 systems in Category 1 of Table 1, and we show 
an example in panel A of Fig. 1. 

(B) Uncertain. For 32 EBs, both eclipses are observed 
but one or more of their measured properties are 
uncertain by more than 20%. This is because one of 
the eclipses is too shallow and/or too narrow. We group 
these 32 systems in Category 2 of Table 1. In panel B of 
Fig. 1, we display an example of a long-period P « 45 day 
EB with a secondary eclipse at $2 ~ 0.55 that is too 
narrow to be accurately measured. 

(C) Roche-lobe filling. Three EBs have wide eclipses 
such that one or both components of the binary must 
be filling their Roche lobes. We list these three systems 
in Category 3 of Table 1, and we show an example in 
panel C of Fig. 1. 

(D) Ambiguous Periods. The orbital periods of 23 of 
our EBs are ambiguous. These 23 EBs can either have 
twin components q ~ 1.0 in nearly circular orbits e ~ 0.0 
or have half the listed orbital periods and exhibit only one 
eclipse per very eccentric orbit. Using the orbital periods 
listed in the OGLE-III LMC catalog, these EBs have 
primary and secondary eclipses that are nearly identical 
and separated by almost precisely 50% in orbital phase. 
Quantitatively, we identify these systems to have values 
of and uncertainties in eclipse depths, widths, and phases 


that satisfy: 


A/-! - A/ 2 | 

< 3[(iJa/i) 2 + (cta/ 2 ) 2 ] ^ 

(2a) 

101-021 

< 3[(cr ei ) 2 + (a 02 ) 2 ] 1/2 

(2b) 

|$ 2 - V2| 

3 (J <J>2 

(2c) 


Given the sensitivity of the data, the observed properties 
imply the 23 systems have large mass ratios q > 0.9 
with extremely small eccentricities e < 0.05 (see §3). 
However, none of the EBs in our sample have eclipse 
depths that satisfy Eqn. 2a (q > 0.9) with secondary 
eclipse phases 3o$ 2 < |$2 — 1 /^\ < 10er$ 2 (e ~ 0.05-0.10). 
Similarly, there is only one EB that satisfies Eqn. 2c (e 
< 0.05) with primary and secondary eclipse depths that 
are discrepant at the (3- 10)cr level (q ~ 0.8- 0.9). Hence, 
there are no twin systems in slightly eccentric orbits, and 
there is only one moderate-mass companion in a nearly 
circular orbit. The prevalence of 23 twin systems in 
nearly circular orbits at these moderate orbital periods 
is therefore highly unlikely. If there is indeed an excess 
of twins in circular orbits relative to twins in eccentric 
orbits, our study does not include them. We expect only 
a few of the 23 EBs that appear to be twins in circular 
orbits to have the listed orbital periods. The majority 
of these EBs more likely have orbital periods that are 
half their listed values, and would therefore exhibit only 
one eclipse per orbit similar to the systems discussed in 
(A) above. In panel D of Fig. 1, we show one example 
where we fold the photometric data with the listed orbital 
period (in black) and the more plausible scenario that the 
binary has half the catalog orbital period (in red). We 
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list these 23 EBs in Category 4 of Table 1. We further 
motivate the removal of these 23 systems in §4 when we 
show the intrinsic frequency of q > 0.6 companions with 
e < 0.2 is relatively sparse. 

(E) Superimposed Variability. Fifteen of the EBs 
exhibit superimposed variability. Three of these systems 
are intrinsic variables, two of which (ID-7651 and 
ID-22929) were already listed as such in the OGLE-III 
LMC EB catalog. The intrinsic variability is readily 
apparent in the unfolded light curves. Moreover, 
the measured intrinsic scatter outside of eclipses is 
substantially higher than the photometric errors, e.g. 
f„j « 2.8 for ID-3414. We note that a few additional 
systems with f a j s=s 1.5-1.9 may exhibit low-amplitude 
variations SI < 0.01 mag, but these variations are 
sufficiently small so as to not to interfere with the 
light curve modeling. We list the three systems that 
exhibit definitive intrinsic variability in Category 5 of 
Table 1. The other 12 EBs exhibit variability in the 
eclipses themselves, only one of which (ID-17017) was 
identified as such in the OGLE-III LMC catalog. For 
these systems, it is possible that more than two bad 
data points occur near the eclipse. More likely, these 
12 EBs display changes in the eclipse depths and/or 
eclipse phases during the seven years of observations. 
Apsidal motion due to tidal and relativistic effects are 
negligible on timescales dt s=s 7 yrs at these wide orbital 
separations. Such evolution in the eclipse parameters 
are most lik ely caused by orb ital motion with a tertiary 
component (IRappaport et all 120131) . We group these 12 
EBs in Category 6 of Table 1, and we display an example 
in panel E of Fig. 1. 

(F) Peculiar. Finally, two EBs have peculiar light 
curves. ID-343 exhibits a pronounced peak in the folded 
light curve at cf> = 0.8 between eclipses. This peak 
may be caused by ellipsoidal modulation in an extremely 
eccentric orbit. ID-4458, which is shown in panel F 
of Fig. 1, displays a sinusoidal variation between two 
eclipses of comparable depth. ID-4458 may contain a 
hot spot and/or disk, and is similar to the green systems 
in the top left corner of Fig. 3 in Paper II. We list these 
two systems in Category 7 of Table 1. 

After removing these 91 systems, our well-defined 
sample contains 130 EBs. We list these 130 systems in 
Category 8 of Table 1. When necessary, we switch the 
primary and secondary eclipses to ensure A/ L > A /2 in 
our well-defined sample. If the epoch of primary eclipse 
minimum t 0 substantially changed from the catalog value 
in order to satisfy this criterion, we place an asterisk next 
to our value of t 0 in Table 1. 

The 130 EBs in our well-defined sample have 
uncertainties in eclipse depths A/i and A /2 and eclipse 
widths ©i and ©2 that are <20% their respective values. 
The uncertainties in the Gaussian analytic fit parameters 
have been used only to determine which EBs have 
detectable and measurable eclipse properties. These 
uncertainties propagate into our Monte Carlo simulations 
when we calculate the fraction of binaries that produce 
detectable eclipses (see §5). The uncertainties in the 
Gaussian analytic fit parameters are not utilized to 
calculate the uncertainties in the physical properties of 
the EBs. Instead, we implement detailed light curve 
models to measure the values of and uncertainties in the 
physical model properties, which we now discuss. 


3. PHYSICAL MODELS (STAGE II) 

3.1. Algorithm 

The physical properties of EBs are routinely 
measu red by fitting detailed mod e ls to the observational 
data dW ilson fc Devinne vl 1 197 lb IP rsa fc Zwitted 120051 : 

iDevor fc Charbonnead 120061 iKallrath fc Milond 12009(1 . 
Normally, spectroscopic radial velocity observations are 
required to dynamically measure the component masses 
M\ and M 2 . In the modern era of wide-field photometric 
surveys, the discovery of EBs is quickl y outpacing 
our a bility to obtain f ollow-up spectra (IDevor et all 
120081: IPrsa et, all l2011al lbh . For large EB samples, 
the physical properties must be inferred based solely 
on the photometric li ght curves. MS const raints 
(IPrsa fc Zwitteil 120051: IK allrath fe Milond [2009) and 
isochrone fitting (IDevor fe Charbonneaul [20 06) have 
helped ascertain EB properties from the photometric 
data. In general, however, these methods lead to large 
systematic uncertainties and/or solutions that are highly 
degenerate. 

In Paper II, we developed a technique that uniquely 
and accurately characterizes the intrinsic physical 
properties of detached EBs with known distances using 
only the photometric data. The distances to EBs in the 
LMC are known. In fact, we have already utilized the 
observed magnitudes (/) and colors (V — I) to select 
EBs with B-type MS primaries. For detached EBs with 
MS primaries, both components are effectively evolving 
along their respective single-star evolutionary sequences. 
The photospheric properties of the stellar components, 
e.g. effective temperatures Tj and T 2 , radii R\ and P 2 , 
and luminosities L\ and L 2 , therefore depend entirely 
on the age r and component masses Mi and M 2 . The 
systematic uncertainties in the evolutionary tracks are 
relatively small, e.g. ~15% uncertainties in the masses 
and ~30% uncertainties in the ages (see Paper II and 
§3.3 for further justification and a full assessment of the 
uncertainties). We can therefore measure the component 
masses Mi and M 2 and ages r of detached EBs with 
known distances based solely on the observed light curve 
features (Fig. 2). 

In our physical models, detached EBs with B-type 
MS primaries can be uniquely described by nine 
independent properties. These nine physical model 
parameters are the orbital period P, epoch of primary 
eclipse minimum t a , primary mass Mi, secondary 
mass M 2 , age r, inclination i, eccentricity e, argument 
of periastron w, and I-band dust extinction Aj. Given 
the age r and component masses Mi and M 2 of the 
binary, we interpolate the radii Pi and P 2 , surface 
gravities gi and g ?, effective temperatures T) and T 2 , 
and luminosities Li and L 2 from pre-MS and MS 
stellar evolutionary tracks with metallic ity Z = 0.008 
(jTognelli et all I 2 TTTTI iBertelli et a l. 2009) . We then use 
the L MC distance modulus p = 18.5 (jPietrzvnski et all 
120131). dust reddening law E(V — I) = 0-7-Ar 
dCardelli et al.lH989l: iFitzpatrickll 1 9991 iNgeow fe Kanbuil 

120051) . and temperature -dependent color indices and 
bolometric corrections ([Pecaut fe Mamaieki l2013f ) to 
transform the intrinsic properties of the binary into 
observed magnitudes and colors. Our physical model 
parameter space (Mi, A/ 2 , r, etc.) of EBs with detached 
configurations, pre-MS/MS evolutionary constraints, 
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Fig. 2.— For detached EBs with known distances and MS 
primaries, the nine observed photometric light curve parameters 
(left) provide unique solutions for the nine independent intrinsic 
physical properties of the system (right). Other properties of the 
binary, e.g. stellar radii R\ and R 2 and luminosities L\ and L 2 , are 
utilized as intermediate indicators (middle), but depend entirely on 
the independent properties Mi, q = M 2 /Mi, and r according to 
stellar evolutionary tracks. 


and known distances is quite different than the typical 
EB parameter space ( T 2 /T 1 , (R 1 +R 2 )/a, etc.) where the 
distances and evolutio nary status o f the compone nts are 
unkno wn (e.g. iDevor fe Charbonneaul 120061 iPrsa et all 
I 2011 af) . 

Using the physical properties of a binary, e.g., P, Mi, 
M 2 , Pi, P 2 , Ti, T 2 , e, ui, etc., we synthesize photometric 
light curves with the EB modeling package NightfaliQ. 
We use the same Nightfall model options adopted in 
Paper II, e.g. a square-root limb darkening law, default 
gravity darkening coefficients, model atmospheres, etc., 
except for three notable distinctions. First, we do not 
assume circular orbits for our EBs at longer orbital 
periods, but instead solve for both the eccentricity e 
and periastron angle w. Second, we set the albedo of 
the secondary to A 2 = 0.7 and implement one iteration 
of reflection effects. Considering reflection effects are 
minuscule for our wider EBs in this study, different 
treatments of reflection have negligible effects on the 
synthesized light curves. Finally, we simulate an EB light 
curve at 1,000 uniformly-spaced discrete orbital phases 
to ensure narrow eclipses are sufficiently sampled. 

Most of our EBs with intermediate orbital periods have 
eccentric orbits and narrow eclipses. Nightfall and all 
other EB software packages that account for tidal effects 
are computationally expensive for eccentric binaries. 
This is because the three-dimensional photospheric 
surfaces of the stars need to be recalculated at each 
of the 1,000 discretely sampled orbital phases. We 
therefore adapt our algorithm from Paper II to guarantee 
fast, automated convergence. Namely, we choose initial 
values for our nine physical model properties that 
are sufficiently close to the true values to ensure y 2 
minimization converges quickly to the global solution. 


2 http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/ 
Nightfall.html 


The major goal of our algorithm is to synthesize light 
curves with Nightfall as few times as possible. Our 
routine can easily be adapted for any population of 
detached EBs with known distances, and can be used in 
combination with any EB light curve modeling software. 

We decompose our algorithm into three steps0 In 
Step 1, we select initial values for our nine physical model 
properties based on the observed light curve features 
quantified in §2. In Step 2, we make small adjustments 
in the physical model properties until the analytic model 
parameters of the synthesized light curve matches those 
of the observed light curve. In Step 3, we utilize a 
Levenberg-Marquardt technique, as done in Paper II, 
to minimize the % 2 statistic between the observed and 
simulated light curves. We elaborate on these three steps 
below. To help illustrate this procedure, we display in 
Fig. 3 the light curve of an example EB, ID-2142, and 
the solutions at the end of each of these three steps. 

Step 1. We use the eight analytic model parameters 
(P, t 0 , (I), Aii, ©i, 4 * 2 , A 1-2, © 2 ) and observed color 
(V — I) from Table 1 to estimate initial solutions for 
the nine physical model properties. In Fig. 2, we show 
how the nine observed light curve features can be used to 
approximate the nine physical properties of the binary. 
We select the physical parameters P and t 0 to match 
the analytic model values. We then estimate e and w 
according to the observed phase of th e secondary eclipse 
and the difference in eclipse widths (IKallrath &; Milonel 
(20091 Eqn. 3.1.24 and 3.1.26, see our Fig. 2): 

7r 

ecosw ~ — ($2 — 72 ) (3a) 

@2 — @1 f ^ 

e since « -——— (3b) 

+ “l 

In this study, u> = 90° if periastron coincides with the 
observed primary eclipse. For our example ID-2142, 
©i « @2 and $2 ~ 0.16, indicating u « 180° and 
e ~ 0.5-0.6. 

The intrinsic colors of B-type MS stars span a narrow 
i nter val —0.3 < (V — I) 0 < —0.1 (IPecaut fe MamaieH 
I20l3ll . We therefore initially assume the intrinsic color 
of an EB to be: 

(V-I) 0 m -0.22 + 0.08((I) - 17) (4) 

where we have accounted for the fact that more luminous 
B-type MS stars tend to be more massive, hotter, and 
bluer. The dust extinction Ai is simply estimated 
from the observed color (V — I) and our adopted dust 
reddening law E(V — I) = (V — I) — (V — I} 0 = 0.7 Aj 
(see Fig. 2). 

We then use the eclipse depths Ap and A / 2 to 
approximate the mass ratio q = M 2 /Mi. For a MS + MS 
binary in a circular orbit, the ratio of eclipse depths 
A / 2 / Aii provides an accurate indicator of the luminosity 
contrast L 2 /L 1 (Fig. 2). This luminosity contrast can 
then be used to infer the mass ratio q according to a 
MS mass-luminosity relation (see Fig. 3 in Paper I). For 

3 The three steps discussed in this section are not to be confused 
with the three full stages of our automated pipeline, which classifies 
EBs (Stage I - §2), fits physical models to the light curves 
(Stage II - §3), and corrects for selection effects (Stage III - §5). The 
three steps regarding physical models are all included in Stage II. 
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eccentric orbits, however, the eclipse depth ratio can be 
modified because the projected distances during primary 
and secondary eclipses can be different. Nonetheless, 
deeper eclipses still suggest larger mass ratios. For 
example, A I 2 > 0.4 mag requires q > 0.7, regardless 
of the eccentricity or whether the secondary is a MS 
or pre-MS star. We use a linear combination of these 
methods to estimate the mass ratio: 


Ain 

q ss O. 6 AI 1 + 0.5A I 2 + 0.5—p- (5) 

Aii 

where the eclipse depths are in magnitudes. 

We next use the observed mean magnitude (I) and 
sum of eclipse widths 0i + 02 to simultaneously measure 
the primary mass Mi and age r. Assuming non-grazing 
eclipses and standard limb darkening coefficients, the 
sum of eclipse widths 0i + ©2 directly provides the 
relative sum of the radii (Ri+R, 2 )/a. Our EBs occupy 
a narrow range of magnitudes 16.0 < (I) < 17.6 
and therefore span a small interval of total masses 
M = Mi + Mi- The orbital separation a oc P 2 / 3 M 1 / 3 
therefore derives mainly from the known period P. We 
can now use 0i + ©2 and P to determine R 1 +R 2 . For 
EBs with B-type MS primaries, we find the following 
approximation: 


R\ + i?2 


0! + q 2 / p y/» 

0 0.01 V 30 days/ 


( 6 ) 


Given our estimates for q and A/ above, we interpolate 
the stellar evolutionary tracks to determine the primary 
mass Mi and age r that reproduce the sum of the radii 
R\ + R -2 according to Eqn. 6 and the observed combined 
magnitude (/). Although (I) and ©i + ©2 both depend 
on Mi and r, they are sufficiently non-degenerate so that 
we can calculate a unique solution. Namely, the primary 
mass Mi largely dictates the luminosity and therefore 
the observed magnitude (I), while the age r primarily 
determines the radii i?i + R 2 and therefore the observed 
eclipse widths ©i + ©2 (see Fig. 2). 

Finally, we select an inclination i that approximately 
reproduces the observed primary eclipse depth A/i 
(Fig. 2). From our estimates of Mi, M 2 = qMi, and 
r, we interpolate the radii Ri and R 2 and effective 
temperatures Xj and X 2 from stellar evolutionary tracks. 
In this step only, we ignore limb darkening and colors of 
the two stars, and instead assume the stars are uniformly 
illuminated grey disks (see Paper I). We assume the 
surface brightnesses of the disks are proportional to 
the stellar temperatures, i.e. the Rayleigh-Jeans law, 
because we are observing at relatively long wavelengths 
in the near-infrared I-band. Using these approximations, 
we calculate the eclipsed area A 0 of the primary at the 
time of primary eclipse t 0 based on the observed primary 
eclipse depth: 


Ah 



A 0 T X 

tt{R 2 iTi + R 2 2 T 2 ) 


(7) 


Given the eclipsed area A a and stellar radii Ri and R 2 , 
we then determine the projected distance d 0 between 
the two stars at t a . The actual ph ysical separatio n at 
primary eclipse is already known via (jKallrath fc Milonel 
120091 Eqn. 3.1.36 evaluated at geometric phase 9 = 0°): 


1 — e 2 

r a = a — - ; - 

1 + e smw 


( 8 ) 


where a derives from our estimates of Mi, M 2 , and 
P according to Kepler’s third law, and e and u> are 
approximated from Eq. 3. Hence, the inclination 
simply derives from cos i = d 0 /r 0 . We limit our 
initial approximation of the inclination to the interval 
i = 86.5° -89.5°. 

We now have initial estimates for the nine physical 
model properties. We emphasize that Eqns. 3-7 are 
simple approximations, and that the true values of e, 
w, A/, Mi, M 2 , t , and i may substantially differ from 
the initial values estimated here. We simply use these 
estimates as initial parameters in our fitting routine in 
order to minimize the number of iterations and accelarate 
convergence toward the final solution (see below and 
§3.2). 

In the top panel of Fig. 3, we compare the I-band 
and V-band light curves of ID-2142 to a simulated 
Nightfall model using the values of the nine physical 
model properties at the end of Step 1. The model 
matches key features of the observed light curve, but 
there are three noticeable differences. First, the 
simulated phase of the secondary eclipse does not match 
the observations; recall that Eqn. 3 is an approximation. 
Second, the simulated color is bluer than the observed 
( V—I ), indicating we underestimated the dust reddening 
Aj in our initial step. Finally, the simulated eclipses 
are slightly deeper than the observed because the more 
accurate Nightfall model accounts for limb darkening 
and color effects. This suggests the actual inclination is 
smaller and/or the mass ratio is slightly different. We 
correct for these visible discrepancies in the following 
step. 

Step 2. Using the nine physical properties from 
Step 1, we synthesize an I-band light curve with 
Nightfall. We then fit the simple analytic model of 
Gaussians (Eqn. 1) as done in §2 to the simulated 
Nightfall light curve. In this manner, we measure 
the analytic parameters of the Nightfall model, e.g. 
Aii 

,mod? ®l,modj 4^2,modi etc. 

We adjust the properties in our physical models 
according to the differences between the simulated and 
observed analytic model parameters. The adjustments 
are motivated as follows. If the modeled eclipse widths 
0 i,mod + © 2 ,mod are wider than the observed ©i + © 2 , 
we select a slightly younger age r (and vice versa). We 
increase the dust extinction Aj if the simulated color 
(V — I)mod is too blue. If the modeled primary eclipse 
A/ymod > Aii is too deep while the modeled secondary 
eclipse A/ 2j mod < AJ 2 matches observations or is too 
shallow, we increase the mass ratio q and decrease the 
inclination i. However, if both simulated eclipses are too 
deep (or both too shallow), we only decrease (increase) 
the inclination i. Finally, we adjust e according to the 
position of and differences in the secondary eclipse phases 
$2 and $ 2 ,mod- In this step, we fix P, t 0 , and ui to 
the values determined in Step 1. Finally, we interpolate 
Mi from the stellar evolutionary tracks based on the 
observed mean magnitude (I) and the revised values for 
r, q, and A/. 

When adjusting our physical model properties, we 
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Fig. 3.— Observed and model light curves for ID-2142. We 
compare the I-band (red) and V-band (blue) OGLE-III LMC data 
to the synthesized I-band (black) and V-band (dotted green) light 
curves at the end of the three steps in our automated procedure. 
We display only the interval —0.15 < (/> < 0.35 that encompass 
the eclipses. Note how the physical model parameters vary only 
slightly between our initial estimate and final solution. 

choose step sizes that scale with the differences between 
the observed and simulated analytic model parameters. 
After making these adjustments, we synthesize another 
I-band light curve with Nightfall. We iterate this step 
until all the analytic model parameters of the simulated 
and observed light curves match within a small tolerance 
level. In the middle panel of Fig. 3, we show our solution 
for ID-2142 at the end of Step 2 after five iterations. 
We therefore required only six Nightfall light curve 
simulations during this middle step. 

Step 3. This final step is essentially the procedure 
outlined in Paper II. We calculate the photometric 
correction factors f a ,j and f„y in both bands. Starting 
with initial model properties determined at the end 
of Step 2 , we utilize a Lev enberg-Marquardt technique 
(MPFIT, IMarkwardtll2009[) to minimize the y 2 statistic 
between the simulated and observed light curves. The 
Levenberg-Marquardt MPFIT algorithm operates by 
independently varying each of the nine physical model 
properties from the previous solution. The routine 
then measures the resulting deviations between the data 
and models, and then calculates a new solution. This 
step therefore requires ten Nightfall simulations per 


iteration. As in Paper II, we simultaneously fit the 
I-band and V-band light curves. We clip up to Af c j + 
Nc.v < 3 data points that exceed 4<r from the best-fit 
model. This results in v = A//+ A/y — A/" Ci j — A/" Ci y — 9 
degrees of freedom. 

In the bottom panel of Fig. 3, we display our 
final solution for ID-2142 after four iterations of the 
Levenberg-Marquardt MPFIT routine. We therefore 
simulated light curves with Nightfall a total of 40 
times in Step 3. The physical model properties changed 
only slightly during this final step. In fact, for ID-2142, 
the variations were all within the uncertainties of the 
physical model parameters. We emphasize that Steps 
1 and 2 were crucial in guarenteeing rapid convergence 
toward the final solution in Step 3. Without them, this 
last step would have required many additional iterations 
or may have converged to a local minimum. 

We utilize this automated procedure for all 130 
detached EBs in our well-defined sample. We present 
our fitted model parameters, physical properties, and 
fit statistics for these systems in Table 2. For MS 
binaries in circular orbits, the deeper primary eclipse 
A/i at time t 0 always corresponds to the smaller, cooler, 
less massive secondary passing in front of the larger, 
hotter, more massive primary. For eccentric orbits, 
however, the situation can be reversed depending on 
the combination of e, w, and i. Indeed, for 18 EBs 
in our well-defined sample, we determined solutions 
such that the less massive component was eclipsed at 
time t 0 . To avoid confusion in nomenclature, we list 
properties in Table 2 according to the primary “p” 
and secondary “s” eclipse features. Namely, M p , 7? p , 
and T p correspond to the component that was eclipsed 
at the epoch of primary eclipse t 0 , and M s , 7? s , and 
T s correspond to the component that was eclipsed at 
the secondary eclipse phase $ 2 - In the text, we refer 
to primary mass Mi = max{M p ,Af s }, secondary mass 
M 2 = min{M p ,M s }, mass ratio q = M 2 /Mi, etc. 

We measure primary masses M-\ = 3.6-13.9Mq, 
which nearly encompasses the full mass range of B-type 
MS stars. We determine mass ratios across the 
interval q = 0.20-1.00, which confirms the OGLE-III 
observations are sensitive to EBs with low-mass 
companions. Our measured dust extinctions cover 
Aj = 0.10-0.58 mag, which is consistent with the range 
of extinctions found in Paper II. Finally, we determine 
ages r = 0.5-190 Myr that span more than two orders 
of magnitude. We further discuss the EB physical 
properties, and their interrelations, in §4. 

Eleven of the 130 EBs have modest fit statistics 
X 2 /^ = 1.10-1.14, i.e. probabilities to exceed x 2 of 
p ss 0.01-0.05 given v « 530 degrees of freedom. Seven 
of these EBs are extremely young with estimated ages 
r < 0.8 Myr (IDs 5153, 7560, 10422,13418,16711, 22691, 
and 22764). The components in these EBs have small 
radii, as demonstrated by their narrow eclipses (Eqn. 6 ), 
and are therefore consistent with the zero-age MS. The 
systematic uncertainties in the stellar evolutionary tracks 
are larger at these young ages, especially considering 
some of the secondaries may still be pre-MS stars (see 
Paper II). Three of the 11 EBs with modest fit statistics 
have primaries at the tip of the MS (IDs 91, 20746, and 
21518), as indicated by their wide eclipses. Again, the 
stellar evolutionary tracks are uncertain at the tip of the 


































MS just prior to the rapid expansion toward the giant 
phase. The one last EB with a poor physical model fit 
(ID-17569) has x 2 /f = 1.11 and p « 0.02. Considering 
our large sample of 130 EBs, we naturally expect 1 - 3 of 
these EBs with modest fit statistics. The remaining 119 
EBs in our well-defined sample have good fit statistics 
0.93 < x 2 / u < 1-09. This is testament that the nine 
independent physical model properties can adequately 
describe detached EBs with known distances and MS 
primaries. 

3.2. Comparison between Initial Estimates 
and Final Solutions 

In the following, we compare the initial estimates for e, 
q, and R 1 +R 2 in Step 1 according to Eqns. 3, 5, and 6 , 
respectively, to the final solutions in Step 3 from fitting 
detailed Nightfall light curve models to the data. We 
can then address the systematic uncertainties in our 
initial estimates and further justify the mapping between 
the basic EB light curve parameters to the physical model 
properties. 

For the 130 EBs in our well-defined sample, we 
compare the initial values of the eccentricities e 
determined from the secondary eclipse phases $2 and 
eclipse widths @1 and ©2 (Eqn. 3) to the final 
Nightfall solutions (top panel of Fig. 4). The initial 
estimates agree quite well with the true final values. The 
rms scatter between the two is only Se = 0.03. This 
validates that Eqn. 3 is more than sufficient for starting 
purposes in our fittng routine. We note that the few 
systems that change by more than Ae > 0.07 between 
Steps 1 and 3 have narrow, poorly sampled eclipses so 
that it is more difficult to precisely measure ©i and 0 2 . 

Similarly, in the bottom panel of Fig. 4, we compare 
the initial estimates of the mass ratios q determined from 
the eclipse depths AR and A / 2 (Eqn. 5) to our final 
values obtained from Nightfall light curve fittings and 
X 2 minimizations. Although the population as a whole 
shows rough agreement between the solutions at the ends 
of Steps 1 and 3, individual systems can substantially 
deviate from the initial estimates. For example, an EB 
with an initial estimate of q ss 0.6 may actually have a 
mass ratio anywhere in the interval q = 0.3-1.0. The 
rms deviation between the initial and final solutions is 
Sq = 0.12, or 5q / q~ 20% the respective values. If we had 
randomly chosen mass ratios q in Step 1 while keeping 
the other initial estimates unchanged, the Nightfall 
light curve solutions in Step 3 would still converge to 
the same final values. We simply find that by adopting 
Eqn. 5 in Step 1 to provide initial estimates for q , the 
number of iterations in Steps 2 and 3 are dramatically 
reduced. 

Finally, we evaluate the discrepancies between the 
values of R 1 +R 2 estimated from the sum of eclipse 
widths ©i +©2 and orbital periods P according to 
Eqn. 6 to the final Nightfall solutions. We 
measure an rms deviation of < 5 (i?i+I? 2 ) = 1.2 Rq, or 
( 5 (i?i+i? 2 ) / (R 1 +R 2 ) ~ 20% the respective values. The 
coefficient in Eqn. 6 should therefore be 7.0 ±1.2 R 0 , 
valid only for EBs with B-type MS primaries. As with 
the mass ratios g, the approximations for R 1 +R 2 based 
on the observed light curve parameters are imprecise but 
sufficiently accurate to provide initial conditions for our 
fitting routine. 




Fig. 4.— Comparison between the initial estimates in Step 1 
of physical model properties based on the observed light curve 
parameters to the final solutions in Step 3 derived from fitting 
Nightfall light curve models. Top panel: the initial eccentricites 
e determined from the phase of the secondary eclipse <l >2 and 
eclipse widths 0j and 0i according to Eqn. 3 correspond quite 
well to the true final values. Bottom panel: the mass ratios q 
estimated from the eclipse depths A/i and A /2 according to Eqn. 5 
are approximate but imprecise indicators of the true mass ratios. 
Nonetheless, the initial estimates are sufficient for starting purposes 
in our fitting routine and, on average, dramatically reduce the 
number of iterations. 

As mentioned in §3.1, R 1 +R 2 is primarily an indictor 
of age r of an EB in our sample rather than the 
component masses Mi and/or M2. For example, at 
age r = 5 Myr, the sum of the stellar radii must 
be contained on the interval R 1 +R 2 ~ 4.4-8.7 R 0 
given any combination of Mi and q > 0.2 that satisfies 
our magnitude limits 16.0 < (I) < 17.6 and —0.25 < 
(V — I) < 0.20 and measured range of dust extinctions 
A 1 ss 0.1-0.5 mag. Meanwhile, at age r = 100 Myr, the 
sum of the radii are systematically larger and confined 
to the interval R 1 +R 2 ~ 5.3-12.3 R 0 given the same 
photometric requirements. Hence, EBs in our sample 
with R 1 +R 2 < 5.3 R 0 must be relatively young while 
those with Ri+R 2 > 8.7 R 0 must be relatively old. We 
initially estimated R 1 +R 2 in Step 1 from the observed 
sum of eclipse widths ©i±©2 according to Eqn. 6 . 
Although not as accurate as the final solutions, Eqn. 6 
provides a model-independent measurement of i?i±i?2- 
The sum of radii R 1 +R 2 estimated from Eqn. 6 is 
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Fig. 5.— The approximate eccentricities e vs. approximate sum 
of stellar radii R 1 +R 2 estimated in Step 1 from the basic observed 
light curve parameters. For an individual system, we indicate 
representative error bars 5e ~ 0.03 and 5(Ri+R,2) ~ 1.2 R©. 
Young EBs with r = 5 Myr that satisfy our photometric selection 
criteria must have R 1 +R 2 = 4.4-8.7 R© (dashed blue), while 
older EBs with r = 100 Myr must have R 1 +R 2 — 5.3-12.3 R© 
(dotted red). The values of R 1 +R 2 , which is an indicator of age 
r, and e are anticorrelated at a statistically significant level. This 
anticorrelation is not caused by selection effects, and is a robust 
and model-independent result. 

therefore a robust and model-independent indicator of 
age r. 

In Fig. 5, we display the eccentricities e measured in 
Step 1 from Eqn. 3 as a function of the approximate 
sum of stellar radii R 1 +R 2 estimated in Step I from 
Eqn. 6. Both sets of parameters are model independent 
and based solely on the observed light curve features. 
According to a Spearman rank correlation test, we 
find that the approximate values of e and R 1 +R 2 are 
anticorrelated (p = —0.18) at a statistically significant 
level (p = 0.04). This suggests that EBs with 
larger components, which are systematically older, favor 
smaller eccentricities. We also compare the 19 EBs with 
approximate R 1 +R 2 < 5.3 R 0 , which must be relatively 
young, to the 27 EBs with R 1 +R 2 > 8.7 R 0 , which must 
be relatively old. According to a K-S test, we find these 
young and old populations of EBs have distributions of 
eccentricities that are discrepant with each other at the 
p = 0.02 significance level. The anticorrelation between 
R 1 +R 2 , which is an indicator of age r, and e is therefore 
statistically significant, robust, and model independent. 
In §4, we further investigate this anticorrelation between 
e and r based on the more accurate final solutions 
obtained from the Nightfall light curve models. 

3.3. Uncertainties 

We now analyze the uncertainties in the final solutions 
of our Nightfall light curve models ( see also Pap er II). 
For each system, we utilize MPFIT dMarkwardtl 120091 ) 
at the end of Step 3 (§3.1) to calculate the measurement 
uncertainties. For all 130 well-defined EBs, the nine 
physical model parameters have unique solutions and 
finite measurement uncertainties. Some of the model 
parameters, however, have solutions that are correlated 
with each other. In addition, uncertainties in the dust 
reddening law, stellar evolutionary tracks, bolometric 
corrections, and Nightfall light curve models can lead 


to large systematic uncertainties in the physical model 
parameters. In the following, we fully investigate the 
measurement uncertainties, parameter correlations, and 
systematic uncertainties in the context of a specific 
example EB, ID-2142. We then determine the median 
total uncertainties of each model parameter (Eqns. 9-17) 
for the entire population of 130 well-defined EBs. 

For our example EB, ID-2142, we explore the physical 
parameter space via a Markov chain Monte Carlo 
(MCMC) technique. Starting with our final solution at 
the end of Step 3, we implement a Metropolis-Hastings 
“random walk” MCMC algorithm to generate and 
select steps in our phase space of nine physical model 
parameters. At each proposed step, we synthesize 
a Nightfall light curve model given the proposed 
nine physical model parameters. The probability p oc 
e~ Ax / 2 of accepting the proposed step is determined 
by evaluating the difference in the y 2 statistic between 
the proposed step and the current solution. Obviously, 
if Ay 2 < 0, the proposed step is always taken. 
If the proposed step is rejected, the step length is 
effectively zero, i.e., the previous solution is counted 
again. We generate proposed steps according to a 
Gaussian distribution with a fixed standard deviation 
for each of the nine physical model parameters. We 
choose the standard deviation in the step sizes so 
that approximately one-third of the proposed steps 
are accepted. We simulate 32,000 proposed steps 
and light curves with Nightfall, which exceeds the 
total number of models generated in §3.1 used to fit 
solutions for all 130 well-defined EBs! It is therefore 
quite computationally expensive to calculate robust 
measurement uncertainties and correlations between 
model parameters for an individual EB with this 
MCMC algorithm. The distribution of the ~12,000 
accepted steps and ~ 20,000 repeated solutions provide 
the nine-dimensional joint probability distribution for 
the physical models. For each of the nine physical 
model parameters, we marginalize across the other eight 
parameters to calculate the one-dimensional probability 
density function. We also compute the two-dimensional 
joint probability distributions for each of the 9 C 2 = 36 
parameter combinations. 

In Fig. 6 , we display the one-dimensional probability 
distributions for the nine physical model parameters 
(diagonal panels) and the two-dimensional joint 
probability distributions for the 36 parameter 
combinations (off-diagonal panels). Although some 
of the parameters are mildly to significantly correlated 
with each other, the measurement uncertainties are finite 
for all nine physical model parameters. The MCMC 
technique confirms the uniqueness and non-degeneracy 
of the physical model solutions. Moreover, the 
measurement uncertainties we determined from the 
robust MCMC algorithm are consistent with the 
measurement uncertainties we evaluated with the 
MPFIT routine. We can therefore rely on the MPFIT 
measurement uncertainties we calculated for all 130 
well-defined EBs. 

The uncertainties in the orbital parameters P and 
t 0 are solely due to the measurement uncertainties 
and dictated by the sensitivity and cadence of the 
OGLE-III LMC observations. The solutions for P and 
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Fig. 6.— The probability density functions (diagonal panels) and joint probability distributions (off-diagonal panels) of the nine physical 
model parameters for ID-2142. We emphasize the confidence intervals/regions in each panel only account for measurement uncertainties. 
We compare the best-fit solutions and lcr measurement uncertainties based on the MPFIT routine (black +’s and green intervals) to the 
68% (red) and 95% (blue) confidence intervals/regions determined from our MCMC technique. For each physical model parameter, we list 
the measurement uncertainty for ID-2142, the systematic uncertainty for ID-2142 (if any), and the median total uncertainty for all 130 
well-defined EBs. These panels demonstrate: (1) the solutions are unique, (2) the uncertainties in P, t 0 , e, and uj are small and primarily 
dictated by the sensitivity and cadence of the OGLE-III LMC observations, and (3) the measurement uncertainties for Mi, M 2 , r, i and Aj 
are correlated, but the systematic uncertainties in the bolometric corrections, stellar evolutionary tracks, and dust reddening law dominate 
the total uncertainties. 


t Q are therefore independent of the other seven model 
parameters (note the fairly circular contours in the first 
and second columns of panels in Fig. 6). For ID-2142, 
we measure la uncertainties of SP « 0.0001 days and 
6t Q « 0.004 days. We find the median lcr uncertainties 
for the entire population of 130 well-defined EBs to be: 

(SP) ft 1.4 x 10~ 5 P ft 0.0004 days (9) 

(St 0 ) ft 0.007 days (10) 

Note that our example ID-2142 has slightly smaller 
uncertainties than average because it is relatively bright 
and its eclipses are well sampled. 

As with P and t Q , the uncertainties in e and ui are 
primarily determined by the sensitivity and cadence of 


the OGLE-III LMC observations. The solutions for e 
and u are therefore independent of the other parameters, 
but are slightly correlated with each other (see last two 
rows in Fig. 6). The eclipses are sufficiently sampled to 
easily break this degeneracy. For ID-2142, we calculate 
a 95% confidence interval of ui = 175° -180°. Note that 
we measured eclipse widths ©i = 0.0047 < 02 = 0.0050, 
also indicating ui < 180° according to the approximations 
in Step 1 (Eqn. 3). Based on the Nightfall light 
curve models, we calculate formal lcr measurement 
uncertainties of <5e meas ~ 0.001 and <5w meas ft 1.4° for 
ID-2142. 

For ID-2142 and some other EBs in our sample, 
the measurement uncertainties <5e meas < 0.005 and 
<5w meas < 1-5° are extremely small. Nightfall treats 
each stellar component as a three-dimensional polyhedral 
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mesh with a finite number of flat surfaces. We 
suspect this finite resolution limits the true sensitivity to 
systematic uncertainties of <5e sys « 0.005 and Suj meas ~ 
1.5°. In any case, the measurement uncertainties de meas 
and dw m eas increase and dominate the total uncertainties 
as the eccentricities e decrease. We measure median total 
uncertainties of Se rs 0.02 and Suj rs 4° for e > 0.5, 
Se ~ 0.03 and Suj rs 10° for e rs 0.3, and Se rs 0.05 and 
Suj rs 20° for e rs 0.1. Obviously, the periastron angle w 
is not defined, and therefore not constrained, if the orbits 
are circular. For the entire population of 130 well-defined 
EBs, we find the following relations adequately describe 
the median total uncertainties: 

(Se) rs 0.06 — 0.07e (11) 

<M rs - (12) 

e 

Solutions for the remaining five parameters Mi, M 2 , 
t, i, and Aj are all correlated with each other (see 
Fig. 6). Moreover, unlike P, t 0 , e, and ui, which have 
relatively symmetric Gaussian errors, the probability 
density functions of Mi, M2, r, i, and Aj are mildly 
to significantly asymmetric. The three parameters M 1 , 
r, and Ai are especially correlated along the observed 
magnitude (I). In other words, solutions with more 
massive primaries M 1 require younger ages r and higher 
extinctions Aj to produce the same observed I-band 
flux. The secondary mass M 2 is also anticorrelated 
with r. Finally, the inclination i mildly depends on the 
three parameters Mi, r, and Aj that are significantly 
correlated with each other. 

Although Mi, M 2 , t, i and Aj are correlated with each 
other, there is sufficient information in the observed light 
curves and our constraints (e.g., distance, evolutionary 
tracks, dust reddening law) to break the degeneracies and 
provide unique solutions (see also §3.1). For example, if 
we were to fix the primary mass at Mi = 11.0 M 0 (i.e., 
the 2.5cr upper limit according to the probability density 
function in Fig. 6), the other parameters would converge 
to M2 = 7.1 Mq, t = 3.6 Myr, i = 86.35°, and Ai = 
0.45 mag with a fit statistic that is Ay 2 = 6.7 larger than 
the best-fit solution. For this larger primary mass, there 
is no combination of r and A/ that can satisfactorally 
reproduce the observed magnitude (I) and color (V — I). 
Similarly, if we were to fix the primary mass at Mi = 
9.2 Mq (i.e., the 2.5cr lower limit), the other parameters 
would converge to M 2 = 6.1 M 0 , r = 13.3 Myr, 
i = 86.08°, and A/ = 0.39 mag with a fit statistic that 
is Ay 2 = 5.8 larger than the best-fit solution. In this 
case, the component masses both decrease by rs12% (to 
maintain the same ratio of eclipse depths All / AI 2 ), and 
so the orbital separation a decreases by 4% according 
to Kepler’s third law. The relative sum of the radii 
(i?i + R 2 )/a, which derives directly from the sum of 
eclipse widths 0i + © 2 , is measured to 1% precision in 
our Nightfall light curve models. If a decreases by 4%, 
then Ri + R 2 must also decrease by rs4%. According 
to the MS stellar evolutionary tracks, if Mi and M 2 
decrease by 12%, then the radii i?i and R 2 decrease by 
9% given the same age r = 7.2 Myr. Hence, the age 
must increase to r = 13.3 Myr so that the sum of radii 
i?i + R 2 only decreases by 4%. If the masses decrease, 


the radii decrease, and the age increases, then the 
temperatures Tf and T 2 both decrease according to the 
stellar evolutionary constraints. However, if Ri and Tf 
both decrease, it is difficult to maintain the same values 
of (I) and (V — I) with only one free extra paramter 
Aj. Hence, there is no combination of M 2 , r, and Aj 
that can satisfactorily reproduce the observed values 
of AI1/AI2, ©1 + © 2 , (I), and (V - I) if Mi = 9.2 M 0 . 
This line of reasoning holds for all EBs in our sample, 
and so the physical model parameters will always have 
unique solutions with finite measurement uncertainties. 
For ID-2142, we measure formal 1 a measurement 
uncertainties of 5Mi, meas rs 0.04Mi Ri 0.5 M 0 , 
SM 2 , meas ~ 0.07M 2 rs 0.4 M 0 , <5r mea s ~ 0.25 t rs 1.8 Myr, 
dime™ ~ 0.1°, and 5A/, meas « 0.01A/ rs 0.006 mag. 
We find similar percentage measurement uncertainties 
in these parameters for the 130 well-defined EBs in our 
sample. 

The systematic uncertainties in Mi, M 2 , t, i and 
Aj can be considerably larger and derive from a 
variety of sources. We first investigate the systematic 
uncertainties in the adopted bolometric corrections. 
Our B-type MS primaries and secondaries span a 
large range of temperatures T rs 10,000-30,000 
K and therefore a broad interval of bolometric 
corrections BC = Mhoi — My rs —3.0—0.3 mag 
dPecaut fc Mamaiekl 1201311 . For the hottest stars in 
our sample with T rs 30,000 K and BC rs —3.0, 
the bolometric corrections are u ncertain by 
SB C rs 0.2 mag, i.e. SBC/BC rs 7% dBertelli et al.ll200ll 
iPecaut fc~M amaick 2013). To propagate this systematic 
uncertainty into our solution for ID-2142, we decrease 
the absolute magnitudes of the bolometric corrections 
by 7% and repeat our fitting routine from §3.1. The 
Nightfall light curve models now converge to a final 
solution of Mi = 9.3 M 0 , M 2 = 6.3 M 0 , t = 7.5 Myr, 
i = 86.35°, and Aj = 0.40 mag. The main effect of 
decreasing \BC\ is to decrease the masses Mi and M 2 . 
This is because more of the flux is radiated in the optical 
and so the component luminosities need to be reduced to 
maintain the same observed magnitude (I) . Fortunately, 
the mass ratio q = M 2 /M 1 is not significantly affected 
by the uncertainties in the bolometric corrections. 
The decrease in masses lead to slightly longer ages 
(to maintain the observed eclipse widths), higher 
inclinations (to maintain the observed eclipse depths), 
and lower extinctions (to maintain the observed color). 
The systematic uncertainties in the physical model 
parameters due to the uncertainties in the bolometric 
corrections are therefore (5My bc = O.llMi = 1.1 M 0 , 
SM 2i bc = O.O 8 M 2 = 0.5 M 0 , (5tbc = 0.04r = 0.3 Myr, 
disc = 0.18°, and JA/^bc = 0.1A/ = 0.04 mag. Because 
the primary mass Mi is mainly dictated by the observed 
(I) and the bolometric corrections, we expect similar 
percentage systematic uncertainties in Mi, M 2 , r, i, and 
Aj for the other EBs in our sample. 

We next propagate the uncertainties in the intrinsic 
colors, observed colors, and dust reddening law. 
The uncertainty in the intrinsic colors of B-type 
MS stars are Ri0.01-0.02 mag ()Pecaut fc MamaieS 
l2013fl . the color calibrations of stars in the OGLE-III 
LMC database are also uncertain by rs 0.01-0.02 
mag (|Udalski et al.1 I2008D . and the coefficient in our 
adopted dust reddening law E(V — I) = 0.70 Aj has 
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a ~1 0 % uncertainty (iCardelli et al.lll989HSchlegel et al.1 
119981 : iFitzpatrickl 119991 iNgeow &: Kanbud 1 200 of) . The 
systematic uncertainty in the dust extinction Aj due 
to dust/color uncertainties is therefore 5-4/,dust/coior = 
max{0.02 mag, 0.1A/}. To confirm this estimate, 
we replace the dust reddening law with E(V — I) 
= 0.63A/ in our models for ID-2142 and repeat our 
fitting routine from §3.1. As expected, we measure 
Aj = 0.48 mag, i.e. the dust extinction increased by 
SAj = 0.1 Aj = 0.04 mag, while the other parameters do 
not vary beyond the measurement uncertainties. 

We finally investigate the uncertainties in the stellar 
evolutionary tracks, including the effects of metallicity 
and rotatio n. We replace the Z = 0.008 tracks from 
iBertelli et all (120091) with the Z = 0.006 non-rotating 
models from IGeorgv et all (I2013T ). We refit ID-2142 and 
measure Mi = 10.5 M 0 , M 2 = 6.8 M 0 , r = 9.0 Myr, 
i = 86.38°, and Aj = 0.45 mag. Hence, the 
systematic uncertainties in the stellar evolutionary 
models, including our ability to interpolate between 
the tracks, dominates the uncertainty in the age 
dTt rac k = 0.26r = 1.9 Myr and inclination ditrack = 0.21°. 
We then replace the evo lutionary t racks with the 
Z = 0.006 tracks from IGeorgv et al.l (|2013f l that are 
rotating on the zero-age MS at v/v CT n = 50% the 
critical break-up velocity. We note that «80% of B-type 
MS stars are rotating at v < 0.5u rr it, ~ 250 km s _1 
(lAbt et al.l 120021 : iLevato &; Grossol l2013lh and our EBs 
with intermediate orbital periods P = 20-50 days may 
have tidally evolved toward slower rotational velocities. 
B-type MS stars initially rotating at v/v cr it = 0.5 have 
equatorial radii that are only ns( 3-4)% larger than their 
pola r radii , bu t MS lifetimes tms that are 20% longer 
(IGeorgv et al.1 120131) . It is therefore the differences 
in the evolutionary tracks of stars with rotation, not 
the distortions in their shapes, that can significantly 
affect our model solutions. We refit ID-2142 with the 
rotating non-synchronized stellar models, and measure 
Mi = 10.5 M 0 , M 2 = 6.9 M 0 , r = 8.6 Myr, i = 86.30°, 
and Aj = 0.45 mag. For ID-2142, the differ ences between 
t he n on-rotating and rotating tracks from IGeorgv et al.l 
( 2013 1 are within the measurement uncertainties. This is 
because the tracks with v/v cr n = 0.5 do not significantly 
deviate from their non-rotating counterparts until the 
ages reach r > 0.8tms the non-rotating MS lifetimes. 
For ID-2142 and the majority of EBs in our sample with 
primary ages r < 0.8tms> the uncertainties due to the 
effects of rotation are negligible. For the few systems that 
are extremely young (r < 1 Myr) or old (t > 0.8tms)) 
we expect slightly larger systematic uncertainties in the 
ages and masses. 

By adding the measurement uncertainties and various 
systematic uncertainties above in quadrature, we 
estimate the total median la uncertainties for the 130 
well-defined EBs to be: 


(8Mi) n 

s 0.15Mi 

(13) 

(SM 2 ) n 

s 0.15M 2 

(14) 

s max{0 

.5 Myr, 0.35r} 

(15) 


(Si) ns 0.4° (16) 

(SAj) ss max{0.03 mag, 0.15A/} (17) 

For these five parameters, the uncertainties are 
dominated by the systematic uncertainties in the 
bolometric corrections, dust reddening law, and 
evolutionary tracks. For older EBs with primary ages 
>80% their MS lifetimes, the total uncertainties in the 
masses (5M\) ns 0.2Mi and (8M 2 ) ns 0.2M 2 and ages (St) 
w 0.45r are slightly larger due to the effects of rotation. 

Some of our EB light curve model solutions can be 
biased due to contamination with a third light source, 
e.g., a tertiary companion or a background/foreground 
object along similar lines of sight. In Papers I and II, 
we estimated that only «10% of our B-type MS EBs in 
the LMC can be contaminated by a third light source 
that is bright enough to significantly contribute to the 
systematic uncertainties. Unlike the previously discussed 
sources of systematic uncertainties that contribute to 
all 130 well-defined EBs, contamination by a third light 
source affects only a small subset of our sample. 

In addition to calculating the uncertainties for the 
nine independent physical model parameters, we also 
estimate the uncertainties in the dependent physical 
properties. The total uncertainties in Mi and M 2 
are «15% but slightly correlated with each other (see 
above). The total median uncertainty in the mass 
ratio is therefore (Sq) = max{0.03, 0A2q}. Because 
the quantity (i?i+i? 2 )/a is precisely constrained from 
the observed eclipse widths, the uncertainties in i?i, 
i? 2 , and a mainly derive from the uncertainties in M\ 
and M 2 according to Kepler’s third law. We measure 
(SRj) ns 0.077?! « O.3R 0 , (6R 2 ) « 0.07i? 2 « 0.2 R 0 , 
and (5a) ns 0.06a ss 6 R 0 . Finally, given the 
«(20-30)% uncertainties in the luminosities (primarily 
due to uncertainties in the bolometric corrections) and 
the s=s7% uncertainties in the radii, the uncertainties 
in the temperatures are ~8% according to the 
Stefan-Boltzmann law. Hence, the total median 
uncertainties are (<5Xi) ns 0.08Xj S3 1,500 K and 
(ST 2 ) ns 0.08 T 2 ns 1,100 K. 

4. EB TRENDS 

In our sample of 130 EBs, several trends and 
correlations exist among the nine physical model 
properties. Most of these trends are caused by 
geometrical and evolutionary selection effects in our 
magnitude-limited sample of EBs. We correct for these 
selection effects in the third stage of our pipeline (§5). 
Two correlations, however, are intrinsic to the population 
of binaries with B-type MS primaries. In this section, we 
first discuss these two empirical relations we uncovered 
from the data, and then we explain the trends that are 
caused by selection effects. 

In Fig. 7, we display the measured I-band dust 
extinctions Aj as a function of age r for the 
130 well-defined EBs. These two parameters are 
anticorrelated (Spearman rank correlation coefficient 
p = —0.34) at a statistically significant level (probability 
of independence p = 8xl0 -5 ). We fit a log-linear trend 
to the total population of 130 EBs (green line in Fig. 7): 
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Fig. 7.— Measured I-band dust extinctions Aj and ages r for 
the 130 EBs in our well-defined sample (black squares). We display 
representative uncertainties for two systems in opposite corners of 
this parameter space. We also display our unbiased subsample of 
98 EBs that is relatively free from photometric selection effects 
(filled blue). The dust extinctions clearly diminish with age, even 
within our unbiased subsample, demonstrating the dust content 
in stellar environments systematically decreases with time. We fit 
a log-linear relation to the total population (green) and unbiased 
subsample (red). The latter is an empirical age-extinction relation 
that can be implemented when modeling other stellar populations. 

Aj,total (mag) = 0.39 - 0.071og( 1] ^ ) (18) 

The slope in the above relation may be biased toward 
negative values due to a photometric selection effect in 
our magnitude-limited sample. Specifically, EBs that 
are intrinsically bluer and more luminous systematically 
contain younger, short-lived, more massive primaries. 
These blue, luminous, younger EBs may therefore require 
larger dust extinctions and reddenings to satisfy our 
photometric selection criteria (and vice versa). In Fig. 8, 
we show the measured absolute magnitudes Mj and 
intrinsic colors (V — I) 0 as a function of dust extinction 
A/ for our 130 well-defined EBs. We also display our 
photometric selection criteria based on the observed 
magnitudes 16.0 < (I) < 17.6 and observed colors 
—0.25 < {V — I) < 0.20 (green lines). Indeed, there are 
several intrinsically red, low-luminosity, older EBs with 
Mi ss —1.2 that are in our sample only because they 
have small dust extinctions Ai ss 0.2 mag. If they were 
to have slightly higher dust extinctions, they would fall 
below our selection limit of (I) = 17.6 (see Fig. 8). 

In §5, we account for our photometric selection criteria 
when analyzing all 130 EBs in our well-defined sample. 
Here, we correct for photometric selection effects by 
further culling our sample according to the intrinsic 
properties of Mj and (V — I) 0 . To obtain an unbiased 
subsample, we can choose EBs across any interval 
of Mi and ( V — I) 0 that also satisfies our selection 
criteria on observed magnitudes and colors. To retain 
most of the sample, we select the regions enclosed by 
-2.63 < Mi < -1.35, -0.341 < (V-I) 0 < -0.115, and 
0.13 < Ai (mag) < 0.45 (red lines in Fig. 8). The 98 
EBs that satisfy these extra selection criteria (filled blue 
systems in Fig. 8) represent an unbiased sample relatively 
free from photometric selection effects. 

Even within this unbiased sample of 98 EBs, the 
intrinsic colors ( V — I) 0 and dust extinctions Ai 



Fig. 8.— Measured absolute magnitudes Mj (top) and intrinsic 
colors (V — I) 0 (bottom) as a function of dust extinction Aj 
for our 130 well-defined EBs (black squares; representative errors 
shown for two systems). We also display the limits on observed 
magnitudes (I) and observed colors (V — I) imposed by our 
photometric selection criteria (green lines). The 98 EBs (filled 
blue) that are enclosed by both red regions are relatively free from 
photometric selection effects. Even within our unbiased sample, 
intrinsically bluer EBs that contain hotter, more massive, younger 
primaries favor larger dust extinctions. 

are still anticorrelated (p = —0.25) at a statistically 
significant level (p = 0.02). As can be observed in 
the bottom panel of Fig. 8, there are relatively few 
intrinsically blue systems (V — I) 0 « —0.30 with small 
dust extinctions Ai ss 0.2 mag. Similarly, there are few 
intrinsically redder EBs (V —1) 0 ~ —0.15 with large dust 
extinctions Ai « 0.4 mag. Intrinsically bluer EBs contain 
hot primaries that are systematically more massive, 
short-lived, and younger. Hence, the anticorrelation 
between age r and dust extinction Ai is real. 

In Fig. 8, we also display r and Ai for the 98 EBs 
(filled blue) in our unbiased subsample. Although not as 
prominent, the ages r and dust extinctions Ai for the 
98 EBs in our unbiased sample are still anticorrelated 
(p = —0.23) at a statistically significant level (p = 0.02). 
For example, there is a complete absence of EBs with 
Ai < 0.2 mag at r < 15 Myr. In contrast, there are 
many EBs in our unbiased sample with A/ = 0.1 - 0.2 mag 
at r > 15 Myr. The intrinsic anticorrelation between 
Ai and r in our unbiased sample demonstrates a 
relationship between dust content and ages of stellar 
environments. Young EBs, and young B-type MS stars 
in general, with r ss 1 Myr are embedded in dusty 
envelopes and/or molecular clouds with photometric 
extinctions A/ ss 0.33 mag. Meanwhile, older EBs with 
r ~ 100 Myr reside in less attenuating environments with 
Ai « 0.26 mag. We fit a log-linear trend to the unbiased 
sample of 98 EBs: 
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Fig. 9. — Distribution of I-band dust extinctions Aj for the total 
population of 130 EBs (black), 42 EBs with ages r < 12 Myr 
(red), and 88 EBs with r > 12 Myr (blue). Although the total 
population peaks at Aj « 0.25 mag with a long tail toward 
high dust extinctions, the young and old subsamples can each be 
accurately described with Gaussian distributions (dotted) centered 
at Aj « 0.34 mag and Aj « 0.25 mag, respectively. 


^/.unbiased (mag) = 0.33 - 0.03 l°g( 1 |^ ) ( 19 ) 

valid for 0.5 Myr < r < 200 Myr (red line in Fig. 8). Even 
after accounting for selection effects, the value of and 
measurement uncertainty in the slope —0.029 ±0.011 
is still inconsistent with zero at the 2.6cr confidence 
level. This is a similar probability of significance based 
on the Spearman rank test above (probability of no 
correlation p = 0.02 between Aj and (V — I) 0 ). The 
p»30% systematic uncertainty in the ages r and «10% 
systematic uncertainty in the extinctions Aj propagate 
into Eqn. 19. The values of and total uncertainties are 
therefore —0.029 ±0.014 for the slope and 0.33 ±0.04 
mag for the mean dust extinction at r ps 1 Myr. The rms 
in the measured dust extinctions Aj around the above 
relation is a = 0.08 mag. 

It had been previously known that younger early-type 
stars in the LMC experience sligh tly highe r dust 
extinctions t h an l ate-type stars dZaritskvl 119991 : 
iZaritskv et al.l I2004H . In the present study, we 
have measured the relationship between age r and 
dust extinction Aj. Quantifying age-dependent dust 
extinctions is crucial when analyzing the spectral energy 
distributi ons of unresolved stellar populations i n distant 
galaxies dPanuzzo et al.ll2007t iSilva et al.ll2011l) . Young 
O- and B-type stars, which dominate the ultraviolet 
component in star-forming galaxies, will experience 
systematically higher dust extinctions than the older, 
redder stars. To accurately constrain the star-formation 
histories of these galaxies, it is imperative to account 
for age-dependent dust extinctions. We note that 
different galaxies and stellar populations will have 
slightly different dust extinctions as a function of age. 
Nonetheless, our empirical age-extinction relation (Eqn. 
19) can provide insight when calibrating models of 
un reso l ved ste llar populations. 

IZaritskv et al.1 (|2004f ) found that the dust extinction 
distribution toward young, hot stars in the LMC peaks 
at Aj ps 0.25 mag with a long tail toward higher 



Fig. 10.— Measured eccentricities e and ages r for the 130 
EBs in our well-defined sample (black squares; representative 
uncertainties shown for two systems). Binaries with B-type 
MS primaries and intermediate orbital periods are preferentially 
born with large eccentricities, which suggest they formed via 
dynamical interactions and/or tidal capture. Moreover, the 
observed slope (red line) in the age-eccentricity anticorrelation 
provides a constraint for dynamical tides in hot MS stars with 
radiative envelopes. 

values. This is consistent with our total population of 
130 EBs with B-type MS primaries (see Fig. 9). By 
dividing our EB population into young (r < 12 Myr) 
and old (r > 12 Myr) subsamples, we find that both 
subsamples can be fitted with simple Gaussians centered 
at Aj ps 0.34 mag and Aj « 0.25 mag, respectively. 
Hence, the non-Gaussian distribution of dust extinction 
may be simply due to a selection effect with age. The 
very young EBs, which represent a small fraction of 
the total population, occupy the long tail toward large 
dust extinctions. Meanwhile, the long-lived EBs, which 
comprise the majority of the sample, form the peak in 
the distribution at Aj ps 0.25 mag. 

We now examine the second physically-genuine trend 
in our EB population. In Fig. 10, we show the measured 
eccentricities e as a function of age r for the 130 EBs 
in our well-defined sample. The eccentricities and ages 
are anticorrelated (Spearman rank correlation coefficient 
p = —0.39) at a statistically significant level (probability 
of no correlation p = 5xl0~ 6 ). 

This observed anticorrelation is primarily because 
eccentricities decrease with time due to tidal evolution. 
The observed trend may be accentuated by a secondary 
effect, whereby EBs with more massive, short-lived 
primaries favor larger eccentricities. However, this 
relation between primary mass M\ and eccentricity e 
cannot fully explain the observed anticorrelation between 
r and e. For example, the eccentricities and ages of the 32 
EBs with massive primaries M\ ps 8.5-13.9 Mq are still 
anticorrelated (p = —0.44) at a statistically significant 
level (p = 0.01). Similarly, the 98 less massive EBs with 
Mj ps 3.6-8.5 M 0 have eccentricities and ages that are 
anticorrelated (p = —0.28) at a statistically significant 
level (p = 0.005). Although EBs with early-B primaries 
may be born with systematically larger eccentricities, 
the anticorrelation between age r and eccentricity e is 
dominated by tidal evolution and is observed in both 
early-B and late-B MS subsamples. 

For late-type stars with M < 1.3 M 0 , orbital energy 
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is most efficiently dissipated into the interior of the 
stars via convective eddies in the stellar atmosphere s 
(fZaEnHT977t Hhffl HMD: IZahnl fl989t IHurlev et, al.ll200l . 
This equilibrium tide model for convective damping 
has been tested against observations of late-type 
bi naries in variou s envir onments with different ages 
(jMeibom fc Mathieull2005f) . For more massive stars M 
> 1.3 M 0 with radiative envelopes, such as our B-type 
MS stars, tides ope rate dynamically via oscillatio ns in 
the stellar interiors (IZahnlfl975l IHurlev et al.ll2002f) . By 
estimating the ages of 130 early-type EBs, we have 
measured the evolution of binary eccentricities due to 
dynamical tides with radiative damping. 

The slope of the observed age-eccentricity 
anticorrelation provides insight into the tidal evolution 
of highly eccentric binaries. We fit a log-linear trend to 
the observations (red line in Fig. 10): 

«,=0.53-0.141og( T5 y (20) 

The value of and measurement uncertainty in the slope 
is —0.14 ±0.03. Hence, the slope is negative at the 50- 
confidence level, similar to the statistical significance 
determined from the Spearman correlation test above. 
Again, systematic uncertainties in the ages r and 
eccentricities e contribute to the uncertainties in the 
coefficients in Eqn. 20. After calculating the total 
uncertainties, we find the mean eccentricity at r « 1 Myr 
is 0.53 ±0.05 while the slope is —0.14 ±0.05. The rms 
scatter in the measured eccentricities around the above 
relation is cr e ss 0.16. 

The intercept in Eqn. 20 implies a circularization 
timescale of r c i rc ~ 5 Gyr for our EBs with B-type 
MS primaries and moderate orbital periods P « 
20-50 days. However, tidal damping is no t as efficient 
when the orbits become less eccentric (tlbul HMDl . 
The true circularization timescale may therefore be 
longer if the age-eccentricity relation flattens beyond 
r > 200 Myr. Conversely, older EBs have systematically 
larger components (see Fig. 5), and so tidal damping 
may become more efficient as the primary fills a larger 
fraction of its Roche lobe. In any case, these short-lived 
B-type MS primaries will expand beyond R\ > 10 R 0 and 
evolve toward the giant branch long before the orbits are 
completely circularized. 

Our young EBs with generally large eccentricities 
experience extreme tidal forces. In fact, a few of the 
EBs with e > 0.6 in our sample have modest Roche-lobe 
fill-factors RLFF ss 0.3 at periastron. Tidal evolution of 
highly eccentric binaries is quite complicated, especially 
considering second-order effect s and non-linear terms can 
become quite important ( Hu tl 119811 b A full analysis of 
tidal evolution in our EB sample is therefore not within 
the scope of the present study. Nonetheless, the observed 
age-eccentricity anticorrelation provides a constraint for 
models of tidal evolution in highly eccentric early-type 
binaries. 

In Fig. 11, we display the cumulative distribution 
function of the eccentricities for the 128 EBs with 
e < 0.68 (green). We do not consider the two EBs 
with e = 0.71 and 0.77 because highly eccentric binaries 
are not complete in our EB sample (see below and 
§5). Moreover, as discussed above, binaries with 



Fig. 11.— Cumulative distributions of eccentricities e for all 128 
EBs with e < 0.68 (green) and subsamples of 91 old EBs with 
t > 10 Myr (red) and 37 young EBs with r < 10 Myr (blue). The 
young population is fully consistent with a thermal eccentricity 
distribution (dashed black), indicating early-type binaries at 
intermediate orbital periods were dynamically captured. 

P = 20-50 days and e = 0.7- 0.8 nearly fill their Roche 
lobes at periastron, and are expected to evolve toward 
smaller eccentricities on rapid timescales. In Fig. 11, 
we also divide our sample into the 91 old EBs with 
t > 10 Myr (red) and 37 young EBs with r < 10 Myr 
(blue). Using a maximum likelihood method, we fit a 
power-law eccentricity probability distribution p e oc e v to 
the observed EBs. We measure 77 = 0.1 ± 0.2, —0.1 ± 0.2, 
and 0.8 ±0.3 for the total, old, and young EB samples, 
respectively. Our total population of EBs (?y = 0.1 ± 0.2) 
is consi s tent w ith the flat distribution (77 = 0 ) observed 
bv lAbtl (I2005T) for his sample of binaries with B-type MS 
primaries and intermediate orbital periods. 

If the orbital velocities and energies of a binary 
population follow a Maxwellian “thermal” probability 
distribution, then the eccentricity probability 
distribution p e = 2e de will be we ighted toward large 
eccentricities (1 Ambartsumian! U937D . Such a population 
of eccentric and thermalized binaries would suggest the 
binaries formed through dynamical interactions, either 
through tidal / disk capture, dynamical perturbations 
in a dense cluster, three-body exc hanges, and/or Kozai 
cycle s with a tertiary companion (iHeggie jl975yPringle 
1986; iTurner et al.1119951 : KroupaUl995t Kiseleva et al. 
1998 : INaoz fc Fabrvckvl 20141) . Surprisingly, the 
observed population of 37 young EBs (77 = 0.8 ±0.3) is 
fully consistent with a thermal eccentricity probability 
distribution (77 = 1; dashed black line in Fig. 11). This 
indicates that massive binaries with intermediate orbital 
periods formed via dynamical interactions on rapid 
timescales r < 5 Myr. 

Previous _ obser vations of spectroscopic 

( Duauennov fc Mavoil 1199 111 and visual 
( Harrington fc Miranianl 119771) solar-type binaries 
have indicated a thermal eccentricity distribution. 
However, these studies recovered the thermal eccentricity 
distribution only after applying large and uncertain 
correction factors for incompleteness. In both the 
spectroscopic and visual binary surveys, the raw 
samples were weighted significantly toward smaller 
eccentricities relative to the thermal distribution. 
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Fig. 12.— Measured eccentricities e versus mass ratios q = 
M 2 /M\ for the 130 well-defined EBs (squares; representative errors 
shown for two systems). There is no evidence for a statistically 
significant correlation between q and e in our sample. We removed 
23 EBs that have ambiguous orbital periods (see item D in §2). If 
we were to fit these 23 systems using the listed periods, they would 
all have q > 0.84 and e < 0.08 (blue diamonds enclosed within red 
lines). Such a dense population in this corner of the parameter 
space is highly unlikely, so it was justifiable to exclude these 23 
systems from our well-defined sample. 

In addition, more recent and complete observ ations 
of solar-type (lAbtl 120061 : IRaehavan et al.l 120101) and 
early-type (jAbtl 120051 ) binaries at intermediate orbital 
periods have revealed a uniform eccentricity distribution 
that is clearly discrepant with a thermal distribution. 

Our raw sample of young early-type EBs is only slightly 
biased toward small eccentricities. In fact, the small 
excess of young EBs with e « 0.1-0.3 relative to the 
thermal distribution in Fig. 11 would be reduced after 
correcting for selection effects. In other words, we expect 
even better agreement between our sample of young 
early-type EBs and the thermal eccentricity distribution 
after considering observational biases (§5). By choosing 
only the EBs with young ages, we have probed the initial 
binary properties of massive stars shortly after their 
formation. For the first time, we have directly observed 
the theoretical thermal eccentricity distribution before 
tides have dramatically reduced the eccentricities. 

In the following, we compare other physical model 
parameters and examine additional trends that could be 
caused by observational biases. We use these observed 
distributions to further justify our selection criteria in 
§2. We also motivate the necessity for incompleteness 
corrections and Monte Carlo simulations, which we 
perform in §5. 

We display the measured eccentricities e as a function 
of mass ratio q in Fig. 12. A Spearman rank test 
reveals no statistically significant correlation (p = 0.25). 
The mass ratios q of early-type binaries are independent 
of their eccentricities e at intermediate orbital periods 
P = 20-50 days. 

In §2 (see item D), we removed 23 EBs with nearly 
identical primary and secondary eclipses separated by 
~50% in orbital phase. We concluded the majority 
of these systems have half their listed orbital periods, 
and therefore exhibit only one eclipse per orbit. If we 
were to fit physical models to these systems assuming 
the listed orbital periods, they would all have q > 0.84 



Fig. 13.— Measured eccentricities e as a function of periastron 
angle ac for the 130 well-defined EBs (squares; representative 
errors shown for two systems). Note that we display the interval 
—270° < a; < 270°, so that some of the systems are repeated. For 
e > 0.4, the concentration of EBs at uj = 0° and uj = —180° = 180° 
as well as deficit at uj = —90° = 270° are due to geometrical 
selection effects. 

and e < 0.08 (blue diamonds within red region of 
Fig. 12). A concentration of 23 EBs in this small corner 
of the parameter space is highly unlikely considering the 
density of systems in the surrounding phase space is 
substantially smaller. We expect only 3 - 5 of the 23 EBs 
to be twins in nearly circular orbits with the listed orbital 
periods. Given only the photometric data, however, we 
cannot easily determine which of the systems are truly 
twins with small eccentricities and which have half the 
listed orbital periods. In §2, we simply excluded all 23 
EBs with ambiguous periods, and we account for the 
incidental removal of the 3-5 genuine systems in §5. 
We emphasize that most of the 23 EBs with ambiguous 
orbital periods have half the listed values, and therefore 
it was appropriate to remove these systems. 

In Fig. 13, we compare the measured eccentricities e 
to the arguments of periastron u for the 130 
well-defined EBs. Assuming random orientations, the 
periastron angle should be uniformly distributed across 
0° < to < 360°. However, the observed systems are not 
evenly concentrated across all u> and e. We notice two 
observational biases in the data, both of which are due 
to geometrical selection effects. 

First, for modest to large eccentricities e > 0.4, the 
EBs cluster near w = 0° and oj = 180°. In fact, the two 
systems with e ~ 0.7-0.8 have uj ~ 0°. For EBs with 
oj = 90° or oj = 270°, one of the eclipses would occur at 
periastron while the other at apastron. The eclipse at 
periaston would be quite narrow according to Kepler’s 
second law, and may be too narrow to be accurately 
measured given the cadence of the OGLE-III data (see 
item B in §2). If the inclination is not sufficiently close to 
edge-on, e.g. i « 87°, then the eclipse at apastron may be 
too shallow to be accurately measured (again, see item B 
in §2). If the inclination was even smaller, e.g. i ss 85°, 
the projected separation at apastron could be large 
enough so there would be no secondary eclipse. These 
systems would exhibit only one eclipse per orbit such as 
those presented in item A of §2. Considering the above, 
it is extremely difficult to observe and measure highly 
eccentric EBs with eclipses that occur near periastron 
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Fig. 14.— Measured mass ratios q = M 2 /Mi as a function 
of age r for the 130 well-defined EBs (squares; representative 
errors shown for two systems). At young ages r < 2 Myr (left 
of blue line), most EBs are in highly eccentric orbits with e ~ 0.6. 
Because of geometrical selection effects, it is difficult to detect these 
EBs, especially if they have small, low-mass companions q < 0.3. 
At older ages r > 12 Myr (right of red line), the primaries are 
systematically larger. The primary eclipse depths A/i, which are 
largely determined by R 2 /R 1 , are therefore shallower. Given the 
sensitivity of the OGLE-111 data, EBs with low-mass companions 
q < 0.3 become undetectable as the primary evolves toward the 
upper MS. 

and apastron. As the eccentricity increases, well-defined 
EBs are only detected as the argument of periastron 
approaches w ss 0° or u> « 180°. 

Second, there is an overabundance of EBs with to « 90° 
relative to those with w « 270°. Quantitatively, there 
are 90 EBs with 0° < w < 180° and only 40 EBs with 
180° < oj < 360°. These two values are discrepant 
at the 4.4cr level according to Poisson statistics. This 
observational bias is due to our definition of the primary 
eclipse minimum t 0 , which determines the reference 
frame for w. Recall the primary eclipse A/i > A /2 at t 0 
must be deeper than the secondary eclipse. If e > 0.2, 
i < 89°, and the primary M\ > M 2 is eclipsed closer to 
apastron, then the eclipse of the most massive luminous 
component Mi may actually coincide with the secondary 
eclipse A I 2 < A/i. Indeed, we found 18 EBs in such a 
configuration whereby Mi is eclipsed at $2 and M 2 is 
eclipsed at t a (see §3 and Table 2). Sixteen of these 18 
EBs have 0° < w < 180°. If we were to define u> according 
to Mi instead of in terms of A/i, then 74 EBs would have 
0° < ui < 180° and 58 EBs would have 180° < w < 360°. 
These two values are now consistent with each other, i.e. 
they only differ at the 1.4cr significance level. 

As indicated above and discussed in §2, we suspect 
the majority of the 48 EBs we removed in items A and 
B of §2 have e > 0.4 and either 20° < w < 160° or 
200° < w < 340°. We test this hypothesis using the 
statistics of the measured systems in our well-defined 
sample. Of the 53 EBs with e > 0.4, 22 have u> < 20°, 
160° < ui < 200°, or u> > 340°. If these 22 systems are 
complete across the specified intervals of u, which total 
80°, and if the intrinsic distribution of periastron angles 
is uniform, then we expect 22 x 360°/80° = 99 EBs with 
e > 0.4. We detected only 53 EBs with e > 0.4, implying 
99 — 53 = 46 EBs did not satisfy our selection criteria. 
These 46 EBs most likley have secondary eclipses that 
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Fig. 15. — Measured primary masses Mi versus mass ratios 
q = M 2 /Mi for the 130 well-defined EBS (squares; representative 
errors shown for two systems). Massive primaries Mj > 11 M@ 
(above blue line) with luminous companions q > 0.6 are too 
bright to be contained in our magnitude-limited sample. Similarly, 
low-mass primaries Mi < 8 Mq (below red line) with q = 0.2-0.3 
companions are either too faint to satisfy our photometric selection 
criteria and/or too old and large to produce detectable eclipses. 

are too narrow, too shallow, or completely absent. This 
prediction of 46 missing EBs nearly matches the 48 EBs 
we removed in items A and B of §2. This consistency 
further demonstrates that geometrical selection effects 
are understood in our sample and the removal of EBs in 
§2 were well-motivated. 

We compare the mass ratios q to the ages r of our 130 
EBs in Fig. 14. There is a lack of extreme mass-ratio 
binaries q < 0.3 at young (r < 2 Myr) and old (r > 12 
Myr) ages. The former is due to geometrical selection 
effects. At extremely young ages, we have shown 
early-type binaries with intermediate orbital periods 
favor large eccentricities. In fact, the median eccentricity 
of the 16 EBs with r < 2 Myr is (e) = 0.59. At these large 
eccentricities, the eclipse that occurs closest to apastron 
will have a larger projected distance, and may therefore 
have a shallower eclipse (see above). Shallow eclipses 
are easily missed given the sensitvity and cadence of the 
OGLE-III observations, especially if the EB contains a 
small, low-mass companion q < 0.3. 

The bias against low-mass companions q < 0.3 at 
older ages is primarily due to an evolutionary selection 
effect. As the primary evolves and expands, the ratio 
of radii R 2 /R 1 decreases and the primary eclipse depth 
All becomes shallower (see Fig. 5 in Paper I). At 
7 w 15 Myr, only companions with q > 0.3 produce 
eclipses A/i > 0.15 mag that are deep enough to be 
detected given the sensitivity of the OGLE-III data. If 
the primary is near the tip of the MS, then q > 0.45 
is required to produce a visible and well-defined eclipse. 
Considering the above, only EBs with ages r ps 2-12 are 
sensitive toward low-mass companions with q ps 0.2-0.3. 

We compare the primary masses Mi to the mass ratios 
q = M 2 /Mi in Fig. 15. There is a clear observational bias 
such that massive primaries Mi = 12-14M© contain 
only small mass ratios q = 0.2-0.5 while late-B MS 
primaries with Mi = 3.6-4.5 Mq include only large 
mass ratios q = 0.6-1.0. This trend is simply due 
to the magnitude limits imposed by our photometric 
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Fig. 16.— Cumulative distribution of mass ratios q = M2/Mi for 
the 13 EBs with M\ = 8-11 Mq, t = 2-12 Myr, and q = 0.22-0.84 
that are relatively free from selection effects (blue solid line). 
Assuming the mass-ratio probability distribution p q oc q 7 can be 
described by a power-law, we display curves for 7 = —3, —2, —1, 
0, and 1 (dotted black). For the 13 unbiased EBs, we measure 
7 = —1.6±0.4 (dashed red), demonstrating binaries with massive 
primaries Mi « 10 Mq and orbital periods P = 20-50 days are 
weighted toward small mass ratios q « 0.2-0.3. 

selection criteria. Massive MS primaries Mi > 12 M 0 
with luminous q > 0.6 MS companions will be brighter 
than our selection limit of (I) = 16.0. Similarly, low-mass 
primaries with Mi < 5 M 0 will be fainter than our 
detection limit of (I) = 17.6 unless there is a bright 
companion q > 0.6 that increases the total luminosity 
of the system. 

The precise mass versus mass-ratio cutoffs in our 
sample also depend on the age of the binary. For 
example, older primaries with Mi ss 5-7M 0 on the 
upper MS will be bright enough (I) < 17.6 to satisfy 
our photometric selection criteria. As stated above, 
EBs with small mass ratios q = 0.2-0.3 produce visible 
well-defined eclipses with Afy > 0.15 mag only when 
the primary is relatively small and young. However, 
young modererate-mass primaries Mi ss 5 - 7 M 0 with 
low-luminosity companions are fainter than our detection 
limit of (/} = 17.6. Hence, our EB sample is sensitive to 
extreme mass ratios q ss = 0 . 2 - 0 . 3 only if Mi > 7M 0 . 
To be conservative, we consider only the primary mass 
interval Mi = 8 -llM 0 to be sensitive to companions 
across the entire interval q = 0 . 2 - 1.0 (distinguished by 
red and blue lines in Fig. 15). 

Considering the above, our EB sample is relatively 
unbiased across the mass-ratio interval q = 0.22-0.84 
(Fig. 12), age interval r = 2-12 Myr (Fig. 14), 
and primary mass interval Mi = 8 -llM 0 (Fig. 15). 
The 13 EBs that are contained in this cube of the 
three-dimensional phase space therefore represent a small 
subsample relatively free from geometrical, evolutionary, 
and photometric selection effects. In Fig. 16, we display 
the cumulative distribution of mass ratios q for these 
13 EBs in our unbiased subsample. Using a maximum 
liklihood technique, we fit a power-law mass-ratio 
probability distribution p q oc q 7 to these 13 EBs. We 
measure 7 = —1.6 ±0.4, demonstrating binaries with 
massive primaries favor extreme mass ratios q ~ 0.2- 0.3 
at intermediate orbital periods P = 20-50 days. 


We emphasize this statistic is based on the small 
unbiased subsample of the 13 EBs, and therefore valid 
only for early-B MS primaries with Mi ss 10 M 0 . The 
median primary mass in our total sample of 130 EBs 
is {Mi) ss 6M 0 . We therefore utilize our total sample 
to derive more accurate statistics as well as probe 
the companion properties of late-B MS stars. In the 
following section, we correct for selection effects so that 
we can make full use of all EBs in our well-defined sample. 

5. CORRECTED BINARY STATISTICS (STAGE III) 

For the final stage of our pipeline, we recover the 
intrinsic binary statistics and distributions by correcting 
for selection effects. As done in Paper II, we first 
determine the probability density functions that describe 
the nine physical model parameters of our EBs (§5.1). 
We then calculate simple estimates for the detection 
efficiencies (§5.2), and then synthesize a large population 
of EBs via a Monte Carlo technique (§5.3). In §5.4, we 
present our results for the intrinsic binary fraction and 
mass-ratio distribution. 

5.1. Probability Density Functions 

We utilize probability density functions similar to those 
in Paper II. For example, we assume random epochs 
of primary eclipse minima t 0 , and that the logarithmic 
orbital periods log P are uniformly distributed across 
P = 20-50 days (i.e., Opik’s law: lAbtl 11983T) . We select 
primary masses Mi =3-30 M 0 and ages r = 0 - 320 Myr 
according to the initial mass function (IMF) and 
star-formation history, respectively, measured for the 
OGLE-III LMC footprint in Paper II. In short, we fitted 
an IMF slope a = —2.4 consistent with the Salpeter value 
and a star-formation history such that the present-day 
star-formation rate is approximately double the rate at 
earlier epochs r = 40-320 Myr. We assume random 
orientations, i.e. cos* = 0-1 and u> = 0°-360° are both 
uniformly distributed across their respective intervals. 

In the present study, we account for the empirical 
age-extinction and age-eccentricity anticorrelations. 
Given an age r, we select dust extinctions Aj according 
to a Gaussian distribution: 

( [Al — Aj unbiasedf 7 ")] \ / 01 \ 

PAi oc exp -^=T- ) (21) 

for 0 < Ai (mag) < 1 and where A/ iUn bi aS ed (t) and 
< 7 Ai = 0.08 mag derive from the fit to the unbiased 
subsample in Eqn. 19. We also choose eccentricities e 
from an age-dependent Gaussian distribution: 

p e oc exp ( — ——( 22 ) 

for 0.0 < e < 0.8 and where eat (r) and a e = 0.16 
derive from Eqn. 20. Finally, we consider the detection 
efficiencies as a continuous function of mass ratio 
q = 0.2-1.0. 

5.2. Simple Estimates 

Before we conduct detailed Monte Carlo simulations, 
we perform simple calculations to estimate the 
probabilities of detecting EBs with P = 20-50 days. 
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For q = 0.8-1.0 companions, the detection efficiencies 
are primarily dictated by two geometrical selection 
effects. First, the orientations must be sufficiently 
close to edge-on. About 90% of our well-defined EBs 
have i > 86.6°, implying the probability of having the 
necessary inclinations to produce observable eclipses is 
Vi = cos (86.6°) = 0.06. 

Second, EBs with longer orbital periods are more 
likely to be missed. Not only do binaries with longer 
periods require larger inclinations i > 87° to produce 
eclipses, but the eclipse widths can also become too 
narrow to be detected given the cadence of the OGLE-III 
LMC observations. We found 73 well-defined EBs with 
P = 20-30 days. Assuming the intrinsic distribution of 
log P is uniform, then we would expect r^ 92 EBs with 
P = 30-50 days. Our well-defined sample includes only 
57 EBs with P = 30- 50 days, suggesting 35 systems were 
missed due to narrow and/or shallow eclipses. Note that 
this is consistent with the 32 EBs we removed in item B 
of §2 with uncertain eclipse parameters. The probability 
that EBs have orbital periods that are sufficiently short is 
therefore Vp = 130 / (130 + 35) ~ 0.8. Considering these 
two factors, the probability of detecting well-defined EBs 
with q « 0.8-1.0 is Vi Vp = 0.06 x 0.8 « 5% (red line in 
Fig. 17). 

For q = 0.2-0.3 companions, we must also consider 
evolutionary and photometric selection effects. As 
discussed in §4, well-defined EBs satisfy our photometric 
selection criteria and are sensitive to low-mass 
companions only if the primaries are relatively young 
with ages r = 2-12 Myr (Fig. 14) and massive with 
Mi « 8-llMg (Fig. 15). Given a typical MS lifetime 
of tms ~ 30 Myr for Mi i=s 8- 11 M 0 primaries, then the 
probability of having the necessary ages r = 2-12 Myr 
is V T « (12 — 2)/30 ~ 0.3. The smallest primary mass in 
our well-defined EB sample is 3.6Mg. Assuming our 
adopted IMF, the probability that an EB contains a 
massive primary Mi > 8 Mgcompared to the probability 
of having any B-type MS primary with Mi > 3.6 M 0 
is Vmi = 0.2. Combining these additional factors, 
then the probability of detecting well-defined EBs with 
q « 0.2-0.3 is V l V P V T V Ml = 0.06x0.8x0.3x0.2 « 
0.3% (blue line in Fig. 17). 

5.3. Monte Carlo Simulations 

We utilize the same technique from Paper II to correct 
for incompleteness across a continuous function of mass 
ratios q. For a given q 1 we select Mi, r, and A/ from their 
respective probability density functions. If the simulated 
binary does not satisfy our photometric selection criteria, 
we generate a new binary. Otherwise, we keep the binary 
and consider its contribution toward the total number 
A/" s i m of simulated binaries. We then select the other 
physical parameters, i.e. t a , P, i, e, and w, from their 
respective probability density functions. 

With the nine physical model parameters for our 
simulated binary, we synthesize an I-band light curve 
with Nightfall. We match the cadence and sensitivity 
of the OGLE-III LMC survey. Specifically, we sample 
the simulated light curve at (A/j) = 470 random epochs 
and add Gaussian noise according to: 


This equation derives from fitting the relation between 
the I-band magnitudes and corrected photometric 
uncertainties for all 221 OGLE-III LMC EBs with 
intermediate orbital periods P = 20-50 days. 

We then fit our analytic model of Gaussians (Eqn. 1) to 
the simulated Nightfall light curve. As in §2, we utilize 
the MPFIT dMarkwardtl 1200911 Levenberg-Marquardt 
routine to measure the values of and uncertainties in 
the eight analytic model parameters, e.g. A/i, <E> 2 , 
etc. To be considered well-defined, we impose the 
same selection criteria adopted in §2. Namely, we 
require the uncertainties in the eclipse depths A/i and 
A /2 and eclipse widths 81 and ©2 to be < 20 % their 
respective values. We also require the fitted orbital 
periods to be unambiguous according to Eqn. 2. If the 
synthesized binary satisfies these selection criteria, we 
consider its contribution toward the total number A/" we ii 
of well-defined EBs. Using a Monte Carlo technique, we 
repeat the above procedure until we simulate 7V" we ii = 50 
well-defined EBs for each value of q. The probability of 
detecting well-defined EBs is simply V = A/weii / -A/sim- 

By creating this mock data set of EBs, we find the 
measurement uncertainties in the eclipse depths and 
widths calculated by MPFIT are themselves uncertain 
by ~20%. Hence, a simulated light curve with Ri 6 cr 
confidence in the analytic light curve parameters may be 
accidentally rejected, while a system with «4 ct detections 
in the eclipse depths and widths may be included 
as a well-defined EB. To account for this systematic 
uncertainty in our selection criteria, we simulate two 
additional sets of EB populations. We first relax our 
criteria and consider EBs as well-defined if the MPFIT 
uncertainties in the eclipse depths A/i and A / 2 and 
eclipse widths ©i and O 2 are <25% their respective 
values. For our final set of simulations, we impose a 
more stringent requirement that the MPFIT relative 
uncertainties are <15%. For each value of q , we therefore 
simulate a total of 3 x A/weii = 150 well-defined EBs. 
In this manner, we have determined the values of and 
uncertainties in V{q). 

In Fig. 17, we display the probabilities V (and their 
uncertainties) of detecting well-defined EBs as a function 
of mass ratio q. As expected, the ability to detect EBs 
with extreme mass ratios q = 0.2-0.3 is substantially 
smaller than the ability to observe EBs with mass ratios 
near unity. At large q > 0.6, the relative uncertainties in 
the probabilities are 5V/V « 11%, which is only slightly 
larger than that expected from Poisson statistics 150 -1 / 2 
= 8 %. Essentially, the majority of EBs with q > 0.6 
have deep and accurately measured eclipse properties, 
and so the precise definition of our selection criteria 
does not significantly affect which EBs are considered 
well-defined. At the smallest mass ratios q = 0 . 2 , 
however, the relative uncertainties more than double to 
8V/V « 25%. All extreme mass-ratio EBs have shallow 
eclipses (see §3) and are close to the detection limit, 
and so the probabilities V of detecting well-defined EBs 
are more uncertain. The probabilities V we calculate 
from our Monte Carlo simulations are consistent with the 
simple estimates derived in §5.2. This demonstrates the 
selection effects are well-understood and the probabilities 
V are reliably measured. 


07 = [1 + io^- 17 0 > /2 ] x 0.0075 mag (23) 
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Fig. 17.— Probability V of detecting well-defined EBs with 
P = 20-50 days as a function of mass ratio q = M 2 /M \. 
The results of our detailed Monte Carlo simulations (black) are 
consistent with our simple estimates (blue and red). In addition 
to orientation effects, evolutionary and photometric selection 
effects in our magnitude-limited sample substantially reduce the 
detection efficiencies for well-defined EBs with extreme mass ratios 
q « 0.2-0.3. 

5.4. Corrected Binary Fraction 

The intrinsic binary statistics are determined by 
weighting each well-defined EB by the inverse of their 
respective probability V(q) of detection as displayed in 
Fig. 17. The total number of B-type MS stars with 
companions q = 0.2-1.0 at P = 20-50 days is simply 
A4omp = «6,500. Given Vb = 96,000 

B-type MS primaries in our photometric sample, then 
T = Vcomp/Ve = 6,500/96,000 « 6.7% of B-type 
MS stars have companions with P = 20-50 days and 
q = 0.2-1.0. The uncertainty in this fraction derives 
from a variety of sources. First, the predicted number 
•Vo. 2 <q<o .3 ~ 2,000 of low-mass companions with q 
= 0.2-0.3 is relatively large but also uncertain. In 
our sample of 130 well-defined EBs, only 8 systems 
have mass ratios q = 0.2-0.3, and so the measurement 
uncertainty from Poisson statistics is 8 -1 / 2 « 35%. The 
systematic uncertainty at q = 0.2-0.3 is ss25% due to 
the uncertainty in the probabilities V of detection (see 
above). The total relative uncertainty in the number of 
companions with q = 0.2-0.3 is therefore «43%, and so 
M).2<g<o.3 = 2,000 ±900. We repeat this calculation for 
the other mass ratio intervals, and find the total relative 
uncertainty in the number of companions is s=s31%, i.e. 
A/"comp = 6,500 ±2,000. Finally, in our Monte Carlo 
simulations, we account for the removal of the 74 systems 
represented in panels A - D of Fig. 1 from our total initial 
sample of 221 EBs. We did not, however, account for the 
17 EBs represented in panels E - F that exhibited variable 
or peculiar eclipse properties. These systems contribute 
a small relative uncertainty of 17/221 = 8%. Hence, the 
total relative uncertainty in the number of companions 
is ~33%. The fraction of B-type MS stars that have 
companions with P = 20-50 days and q = 0.2-1.0 is 
therefore T = (6.7 ± 2.2)%. 

Surveys for double-lined spectroscopic binaries (SB2s) 
with early-type primaries are generally complete for 
modest mass r atios q > 0.25 and short orbital pe riods 
P < 20 days (lAbt et al.l 119901 iSana et all I2012D . In 



Fig. 18.— Cumulative distribution of mass ratios 
M 2 /Mi = 0.2-1.0 for the 130 EBs in our well-defined 
sample (blue) after weighting each system by the inverse of their 
respective detection probability V{q). Assuming the mass-ratio 
probability distribution can be described by a power-law p q oc g 7 , 
we display curves for exponents 7 = —3, —2, —1, 0, and 1 (dotted 
from top to bottom). After correcting for selection effects, we 
measure 7 = —1.1 ±0.3 (dashed red), demonstrating early-type 
binaries with intermediate orbital periods are weighted toward 
extreme mass ratios. 

a sample of 109 B-type MS stars, lAbt et ahl (|1990T ) 
found seven SB2s with q > 0.25 and P = 2-20 days. 
Similarly, in a sample of 71 O-type stars, [Sana et all 
• 2012 ; identified 18 SB2s across the same mass-ratio and 
period intervals. These statistics imply /i og p = 7 / 109 = 
0.06 ± 0.02 and /i og p = 18 / 71 = 0.25 ± 0.06 companions 
with q > 0.25 per decade of orbital period at logP (days) 
= 0.8 for B-type and O-type stars, respe c tivel y. As 
discussed in Paper II and in iChini etahl ( 2012 ). the 
close binary fraction dramatically increases with primary 
mass. 

Based on our B-type MS EBs, we measure 
a (5.6 ±1.4)% corrected binary fraction across 
P = 20-50 days and q > 0.25. This results in 
/ logP = (0.056 ±0.014)/(log 50-log 20) = 0.14 ±0.04 
companions with q > 0.25 per decade of orbital period 
centered at logP (days) = 1.5. This value is consistent 
with the early-type spectroscopic binary fraction 
measured at short orbital periods P < 20 days, implying 
the intrinsic period distribution of early-type binaries 
closely resembles Opik’s law (!Abtlll983t lAbt et ahlll990H . 

5.5. Corrected Mass-ratio Distribution 

In Fig. 18, we display the cumulative distribution 
of mass ratios q after weighting each well-defined EB 
by the inverse of their respective probability V of 
producing observable eclipses. By fitting a power-law 
probability distribution p q oc q 7 to the corrected 
mass-ratio distribution, we measure 7 = —1.1 ±0.3. 
This is consistent with our estimate in §4 (Fig. 16) of 
7 = —1.4 ± 0.3 based on a relatively unbiased subsample 
of 13 young EBs with early-B MS primaries Mi ss 10 M 0 . 
In both cases, binaries with B-type MS primaries and 
orbital periods P = 20-50 days favor small mass ratios 
q = 0.2-0.3. 

Observations of early-type spectroscopic and eclipsing 
binaries with short orbital periods P = 2-20 days reveal 
a mass-ratio probability distribution that is only slightly 
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weighted toward small values, e.g . 7 ~ —0.9-— 0.2 
(|Abt et all 119901: iSana et all 120121 Papers I and II). 
In addition to the power-law component, close massive 
binaries with P < 2 0 days exhib i t a sm all excess o f twins 
with q > 0.9 (ITokovininl 120001 IPinsonneault fe Stanekl 
I2006L Paper I). The preponderance of close binaries 
with moderate mass ratios and excess of twins suggest 
early-type binaries with P < 20 days coevolved 
via fragmentation and competitive accretion i n the 
circumbinary disk ( Abt et al.l Il9 90l: ITokovininl 120001 : 
iBonnell fc Bate! 1200.4 iKouwenhoven et al.ll2009f) . 

In contrast to close massive binaries with P < 20 days, 
we find early-type binaries at moderate orbital periods 
P = 20-50 days are even further weighted toward 
extreme mass ratios, i.e. 7 = —1.1 ±0.3 for our 
total sample and 7 = —1.6 ±0.4 for early-B MS 
primaries. In addition, there is no evidence for an excess 
population of twins at intermediate orbital periods. 
Previous spectroscopic surveys have indicated that the 
mass-ratio probability distribution becomes weighted 
toward smaller values with increasing 1 orbital period 
(|Abt et all 1 1 9901 [kobulnickv et all 1 2 0141) . However, this 
result is primarily based on the smaller frequency of 
SB2s at intermediate orbital periods, especially when 
compared to the frequency of single-lined spectroscopic 
binaries. In the present study, we have measured 
the mass-ratio probability distribution at intermediate 
orbital periods. Our results indicate that early-type 
binaries at slightly longer orbital periods P = 20-50 days 
have experienced less coevolution. In a future paper 
(Moe et al., in prep.), we will analyze early-type binaries 
discovered through other observational techniques, e.g. 
long-baseline interferometry, adaptive optics, common 
proper motion, etc., and investigate this anticorrelation 
between P and q in a more thorough and self-consistent 
manner. 

6 . SUMMARY 

Eclipsing Binary Sample (§2). We analyzed the 221 
EBs in the OGLE-III LMC database with B-type MS 
primaries and orbital periods P = 20-50 days. After 
fitting analytic models of Gaussians to the observed 
light curves, we identified 130 detached EBs that exhibit 
two well-defined eclipses per orbit. The remaining 90 
EBs have uncertain, peculiar, and/or variable eclipse 
properties, including 12 systems that displayed changes 
in the secondary eclipse parameters most likely due to 
orbital motion with a tertiary companion. 

Physical Models /§ 3). We developed an automated 
procedure to robustly and quickly fit detailed physical 
models to the EB light curves. Our algorithm can 
be adapted for any population of detached EBs with 
known distances and MS primaries. We implemented 
our procedure on our 130 detached well-defined EBs 
to measure their intrinsic physical properties, including 
their ages r, component masses Mi and M 2 , dust 
extinctions A/, and eccentricities e. We incorporated 
various techniques to demonstrate the uniqueness and 
robustness of the model solutions as well as the accuracy 
of the model parameters. 

Age-Extinction Anticorrelation /§/). Even after 

considering selection effects, we find the ages r and 


dust extinctions A/ are anticorrelated (p = —0.23) at 
a statistically significant level (p = 0.02). This suggests 
young stars with r < 10 Myr are embedded in dusty 
envelopes and/or molecular clouds with A/ « 0.35 mag, 
while older stars with r > 100 Myr reside in less 
attenuating environments with Aj « 0.25 mag. This 
empirical relation between r and A/ should prove 
beneficial when modeling stellar populations. 

Age-Eccentricity Anticorrelation /§/). We also 
discover the ages t and eccentricities e are anticorrelated 
{p = —0.39) at a statistically significant level 

(p = 5 x 10” 6 ) due to tidal evolution. The slope 
in the observed trend provides a diagnostic for the 
radiative damping constant via dynamical tides in highly 
eccentric binaries with hot MS components. We note the 
tidal circularization timescales e/e in highly eccentric 
binaries with e ~ 0.5-0.8 may be orders of magnitude 
shorter than the circularization timescales when the 
eccentricities e < 0.4 are already small. 

Initial Eccentricity Distribution /§//. We find that 
massive binaries at P = 20-50 days are initially born 
with larger eccentricities (e) « 0.6. Assuming a 
power-law eccentricity probability distribution p e oc e v , 
we measure 77 = 0.8 ± 0.3 for our young early-type 
EBs with r < 10 Myr. This is consistent with a 
Maxwellian “thermal” eccentricity distribution (77 = 1), 
which indicates massive binaries with intermediate 
orbital periods formed via dynamical interactions, either 
through tidal / disk capture, dynamical perturbations 
in a dense cluster, three-body exchanges, and/or Kozai 
cycles with a tertiary companion. 

Binary Fraction /§5/. After utilizing a Monte 
Carlo technique to correct for selection effects, we 
measure that (6.7 ± 2.2)% of B-type MS stars have 
companions with P = 20-50 days and q = 0.2-1.0. 
The frequency of companions per decade of orbital 
period at log P(days) = 1.5 is consistent with 
spectroscopic observations of close massive binaries at 
log P(days) = 0.8. This suggests the intrinsic period 
distribution of binary companions to B-type MS stars 
closely resembles Opik’s law for P < 50 days. 

Mass-ratio Distribution ($5). In our corrected binary 
sample with B-type MS primaries (Mi) = 6 Mg, we 
measure a mass-ratio probability distribution p q oc q 7 
weighted toward small values (7 = —1.1 ±0.3). There 
is a slight indication that binaries with early-B MS 
primaries (Mi) = 10 Mg are even further skewed toward 
extreme mass ratios (7 = —1.6±0.4). Close massive 
binaries with P < 20 days favor moderate mass ratios and 
exhibit a small excess of twin components q > 0.9. This 
indicates our early-type MS binaries with intermediate 
orbital periods P = 20-50 days have experienced 
substantially less coevolution via fragmentation and 
competitive accretion in the circumbinary disk. 
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Table 1: Analytic model parameters that describe the basic light curve features for the 221 EBs with OGLE-III LMC 
catalog properties 16.0 < (I) < 17.6, —0.25 < (V — I) < 0.20, and P = 20-50 days. Based on the measured analytic model 
parameters, we divide the total sample into eight categories: (1) EBs without secondary eclipses, (2) EBs with uncertain eclipse 
parameters, (3) Roche-lobe filling EBs, (4) EBs with ambiguous orbital periods, (5) intrinsic variables, (6) EBs with variable 
ecl ipses, (7) pec uliar EBs, and (8) detached well-defined EBs. For each category, we list the OGLE-III LMC catalog properties 
dGraczvk et all 1201 il l including the identification number, mean color (V — I), and number Mi of I-band measurements. We 
then list the eight best-fit analytic model parameters: orbital period P (days), epoch of primary eclipse minimum t 0 (Julian 
date — 2450000), mean magnitude (/), primary and secondary eclipse depths AJi and A /2 (mag), orbital phase of secondary 
eclipse $ 2 , and eclipse widths ©1 and ©2 (fraction of the orbital period). Finally, we list the fit statistics, including the correction 
factor J a .i in the photometric errors, number M c of clipped data points, and goodness-of-fit statistic Xg/ 11 - 

Category 1: list of 16 EBs without visible secondary eclipses. These EBs most likely have a certain combination of e, u>, and i 
so there is only one eclipse per orbit. 


Catalog Properties 

Analytic Model Parameters 

Fit Statistics 

ID 

(V-I) 

Mi 

p 

to 

(I) 

Ah 

AI 2 

$2 

0 i 

©2 

f*,i 

Me 

Xc,M 

316 

0.01 

211 

24.9947 

3538.528 

17.34 

0.22 

~ 

- 

0.0074 

- 

1.00 

0 

1.01 

2164 

-0.04 

440 

43.0029 

3569.803 

17.34 

0.40 

- 

- 

0.0029 

- 

1.09 

1 

1.01 

3091 

-0.08 

424 

27.4499 

3602.602 

17.14 

0.21 

- 

- 

0.0054 

- 

1.14 

0 

1.09 

4652 

0.02 

445 

23.7441 

3572.959 

16.61 

0.10 

- 

- 

0.0040 

- 

1.06 

2 

1.03 

5548 

-0.13 

466 

23.5235 

3585.619 

17.23 

0.33 

- 

- 

0.0027 

- 

1.09 

0 

1.02 

5704 

-0.12 

421 

25.2059 

3566.122 

16.25 

0.18 

- 

- 

0.0044 

- 

1.20 

0 

1.04 

5973 

-0.05 

468 

26.6085 

3635.421 

16.29 

0.42 

- 

- 

0.0039 

- 

1.20 

0 

1.08 

9850 

0.03 

422 

37.3082 

3634.111 

17.39 

0.49 

- 

- 

0.0025 

- 

1.11 

1 

1.03 

12084 

-0.03 

477 

29.8000 

3588.434 

16.90 

0.18 

- 

- 

0.0037 

- 

1.15 

0 

1.02 

12913 

-0.08 

493 

25.2768 

3589.836 

16.08 

0.09 

- 

- 

0.0050 

- 

1.23 

0 

1.04 

13991 

-0.02 

449 

30.6988 

3592.784 

17.16 

0.41 

- 

- 

0.0029 

- 

1.08 

2 

1.02 

17232 

-0.03 

454 

26.9222 

3639.974 

16.24 

0.12 

- 

- 

0.0023 

- 

1.37 

2 

1.07 

21007 

-0.04 

367 

33.7279 

3600.454 

16.40 

0.09 

- 

- 

0.0097 

- 

1.14 

1 

1.04 

22467 

-0.14 

437 

20.2896 

3593.730 

16.46 

0.14 

- 

- 

0.0028 

- 

1.06 

0 

1.03 

23086 

0.08 

434 

25.0090 

3608.833 

17.30 

0.21 

- 

- 

0.0017 

- 

1.09 

0 

1.02 

25112 

-0.04 

391 

39.9125 

3581.880 

16.07 

0.09 

- 

- 

0.0083 

- 

1.49 

0 

1.04 


Category 2: list of 32 EBs with uncertain eclipse parameters. These EBs generally have eclipses that are too narrow (0 < 0.003) 
and/or too shallow (A I < 0.15 mag) to be accurately measured. 


Catalog Properties 


Analytic 

Model Parameters 



Fit Statistics 

ID 

(V-I) 

Mi 

P 

to 

(I) 

Ah 

AI 2 

$2 

©1 

©2 

fey,I 

Me 


219 

0.01 

405 

49.9100 

3626.670 

17.53 

0.42 

0.56 

0.390 

0.0090 

0.0011 

1.00 

1 

0.99 

1450 

0.11 

465 

28.5096 

3573.275 

17.52 

0.38 

0.17 

0.426 

0.0045 

0.0048 

1.17 

0 

1.01 

1924 

0.06 

464 

31.0128 

3607.649 

16.89 

0.25 

0.09 

0.651 

0.0053 

0.0070 

1.15 

0 

1.00 

2539 

-0.06 

457 

38.1274 

3585.424 

17.07 

0.14 

0.13 

0.756 

0.0033 

0.0021 

1.03 

1 

1.02 

2843 

-0.05 

440 

26.0647 

3599.963 

17.57 

0.17 

0.46 

0.746 

0.0096 

0.0013 

1.07 

0 

0.98 

3492 

0.06 

457 

20.8220 

3598.229 

17.58 

0.48 

0.09 

0.769 

0.0033 

0.0074 

1.15 

0 

1.02 

3745 

0.03 

225 

26.5908 

3575.121 

17.58 

0.41 

0.53 

0.527 

0.0094 

0.0015 

1.11 

0 

1.08 

4095 

-0.05 

432 

33.9572 

3615.866 

17.39 

0.19 

0.15 

0.936 

0.0061 

0.0022 

1.14 

2 

1.00 

4396 

-0.03 

358 

25.3069 

3598.022 

17.14 

0.14 

0.09 

0.542 

0.0060 

0.0059 

1.20 

1 

1.10 

5257 

-0.02 

463 

35.1506 

3575.663 

17.29 

0.16 

0.07 

0.472 

0.0044 

0.0036 

1.05 

2 

0.99 

6494 

-0.12 

421 

45.3346 

3623.790 

16.96 

0.18 

0.14 

0.606 

0.0030 

0.0017 

1.25 

0 

1.01 

7832 

-0.07 

476 

26.6117 

3613.760 

17.35 

0.12 

0.06 

0.468 

0.0059 

0.0016 

1.05 

0 

1.02 

7954 

-0.05 

435 

27.1161 

3606.847 

17.25 

0.20 

0.32 

0.709 

0.0044 

0.0030 

1.11 

1 

1.00 

8824 

0.01 

437 

37.6744 

3617.323 

17.56 

0.18 

0.18 

0.277 

0.0052 

0.0037 

1.14 

1 

1.06 

10248 

-0.14 

557 

43.9229 

3633.148 

16.36 

0.35 

0.72 

0.258 

0.0047 

0.0005 

1.49 

2 

1.00 

11655 

-0.02 

477 

30.0005 

3595.741 

16.29 

0.10 

0.07 

0.431 

0.0056 

0.0150 

1.23 

0 

1.00 

12065 

-0.08 

459 

29.0408 

3563.836 

16.81 

0.20 

0.21 

0.806 

0.0048 

0.0025 

1.09 

1 

1.06 

12202 

-0.17 

477 

41.4788 

3572.545 

16.76 

0.35 

0.45 

0.217 

0.0048 

0.0013 

1.11 

2 

1.03 

12696 

-0.09 

457 

40.3637 

3587.966 

16.76 

0.14 

0.05 

0.373 

0.0050 

0.0033 

1.46 

1 

1.02 

13076 

-0.08 

493 

33.6509 

3637.571 

16.01 

0.06 

0.03 

0.436 

0.0036 

0.0023 

1.07 

1 

1.01 

14307 

0.01 

540 

22.3483 

3601.494 

17.58 

0.12 

0.23 

0.290 

0.0055 

0.0041 

1.00 

1 

1.01 

16651 

-0.10 

449 

45.3897 

3684.769 

17.54 

0.34 

0.69 

0.549 

0.0045 

0.0015 

1.06 

0 

1.02 

16922 

0.02 

580 

41.6338 

3670.541 

17.58 

0.13 

0.30 

0.186 

0.0037 

0.0018 

1.22 

2 

1.00 

17204 

-0.01 

473 

31.9060 

3650.939 

17.60 

0.23 

0.07 

0.566 

0.0032 

0.0028 

1.06 

0 

1.00 

17262 

-0.14 

626 

33.8546 

3650.508 

16.50 

0.34 

0.19 

0.499 

0.0014 

0.0016 

1.30 

2 

1.02 

17957 

-0.11 

626 

47.9194 

3623.161 

17.23 

0.19 

0.05 

0.492 

0.0030 

0.0050 

1.15 

2 

1.07 

18800 

0.03 

470 

22.0031 

3582.822 

17.47 

0.26 

0.17 

0.325 

0.0044 

0.0050 

1.10 

0 

1.02 

20667 

0.10 

417 

35.3260 

3627.664 

17.54 

0.36 

0.60 

0.498 

0.0022 

0.0014 

1.05 

1 

0.97 

22464 

-0.08 

437 

20.7824 

3575.419 

17.38 

0.30 

0.10 

0.499 

0.0034 

0.0036 

1.00 

2 

0.99 

22512 

0.08 

437 

25.1642 

3588.758 

17.41 

0.22 

0.25 

0.484 

0.0091 

0.0440 

1.31 

0 

0.99 

22853 

-0.10 

430 

23.4011 

3637.436 

16.94 

0.19 

0.06 

0.398 

0.0020 

0.0049 

1.12 

1 

1.02 

23330 

-0.02 

414 

43.4464 

3616.771 

16.94 

0.26 

0.35 

0.529 

0.0035 

0.0029 

1.28 

0 

1.01 
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Table 1 (cont.): 

Category 3: list of 3 Roche-lobe filling EBs, as demonstrated by their wide eclipses 0 > 0.06. 


| Catalog Properties 

Analytic Model Parameters 

Fit Statistics [ 

ID 

(V-I) 


p 

to 

(I) 

A/i 

1 A/2 

<1?2 

©1 

©2 

fa.I 

KSi 

1 Xc> 1 

3864 

-0.09 


26.6937 

3602.081 

16.41 



0.497 

0.0640 

0.0670 

■WK1 

LjJ 


16199 

-0.04 


46.3408 

3632.722 

17.18 



0.496 

0.0609 

0.0668 




25591 

0.18 


20.7526 

3588.389 

16.01 


| 0.31 

0.503 

0.0777 

0.0856 


2 



Category 4: list of 23 EBs with ambiguous orbital periods. These systems have A/i « A/ 2 , $2 « 0.5, and ©1 « ©2 given the 
listed orbital periods. The majority of these EBs most likely have half the listed orbital periods, and therefore exhibit only one 
eclipse per orbit such as the systems listed in Category 1. 


Catalog Properties 

Analytic Model Parameters 

Fit Statistics 

ID 

(v-i) 


p 

to 

(I) 

Ah 

A I 2 


©1 

©2 

fa,I 



675 

0.08 


20.5259 

3584.169 

17.29 

BI1EEB 

0.47 

0.500 

0.0079 

0.0071 



1.05 

885 

-0.04 


23.5967 

3564.996 

17.31 


0.33 

0.500 

0.0028 

0.0027 



1.01 

1478 

0.14 


32.0021 

3628.807 

17.09 


0.14 

0.501 

0.0020 

0.0026 


1 

1.03 

8321 

0.00 


34.8297 

3643.253 

17.12 


0.11 

0.502 

0.0031 

0.0025 


2 

1.02 

8376 

0.01 


36.9128 

3638.843 

17.33 


0.23 

0.500 

0.0025 

0.0029 



1.00 

9146 

-0.04 


21.1796 

3606.499 

17.40 


0.25 

0.500 

0.0023 

0.0021 



1.00 

11930 

-0.24 


42.0006 

3542.808 

16.06 


■IMIM 

0.500 

0.0012 

0.0013 


1 

1.03 

11931 

0.06 


24.6953 

3569.815 

16.67 


0.29 

0.500 

0.0050 

0.0048 



1.05 

14753 

-0.01 


22.4059 

3580.305 

17.41 



0.500 

0.0023 

0.0024 


1 

1.02 

15309 

0.06 


44.9548 

3621.008 

16.54 



0.501 

0.0029 

0.0019 



1.00 

17257 

-0.03 


30.1341 

3586.993 

16.48 



0.500 

0.0016 

0.0019 


mm 

1.00 

17407 

-0.06 


25.8992 

3612.024 

16.92 


0.12 

0.499 

0.0053 

0.0048 


mm 

1.08 

17715 

-0.07 


35.5355 

3616.239 

16.73 


0.13 

0.499 

0.0022 

0.0024 


1 

1.03 

18138 

0.12 


42.9000 

3599.789 

16.92 



0.500 

0.0024 

0.0024 



1.03 

19309 

-0.09 


32.1801 

3593.017 

16.13 



0.500 

0.0015 

0.0018 



1.01 

19582 

-0.11 


20.0314 

3608.284 

16.59 



0.498 

0.0024 

0.0026 


1 

1.03 

19612 

-0.01 


20.8496 

3616.468 

17.14 


0.41 

0.500 

0.0026 

0.0029 



1.03 

19651 

-0.06 


26.4068 

3598.679 

16.83 


0.19 

0.499 

0.0026 

0.0034 


1 

0.99 

20441 

0.08 


23.5402 

3608.235 

17.21 


0.13 

0.498 

0.0042 

0.0051 



0.98 

20661 

0.07 


20.5829 

3618.710 

17.41 


0.13 

0.499 

0.0036 

0.0045 



0.99 

21273 

0.04 


20.7646 

3598.007 

16.77 


0.25 

0.499 

0.0031 

0.0037 



1.02 

21477 

-0.02 


23.3476 

3602.139 

17.52 



0.499 

0.0044 

0.0035 


2 

0.97 

24604 

-0.02 


38.0192 

3591.216 

17.04 


0.17 

0.500 

0.0028 

0.0033 



0.99 


Category 5: list of 3 EBs that are intrinsic variables, as indicated by their large rms scatter f a ,i > 1.6. 


| Catalog Properties 

Analytic Model Parameters 

Fit Statistics [ 

ID 

(V-I) 


p 

to 

(I) 

A h 

1 a/ 2 

<|>2 

©1 

©2 

fa.I 



3414 

-0.12 


38.9097 

3607.475 

17.02 


cngi 

0.897 

0.0041 

0.0021 



1.02 

7651 

-0.01 


34.8819 

3595.802 

16.64 


0.12 

0.267 

0.0075 

0.0038 



0.99 

22929 

-0.06 


26.4805 

3597.217 

16.63 


0.22 

0.772 

0.0068 

0.0064 


2 

1.06 


Category 6: list of 12 EBs with variable eclipses. These EBs either have more than Af c > 2 bad data points near the eclipses 
or, more likely, exhibit variations in the eclipse parameters due to orbital motion with a tertiary companion. 


Catalog Properties 

Analytic Model Parameters 

Fit Statistics 

ID 

(V-I) 


p 

to 

(I) 

A/i 

AI 2 


©1 

©2 

fa,I 

97M 


3112 

-0.12 


20.8143 

3569.074 

16.96 


0.33 

0.498 

0.0023 

0.0012 



1.16 

3233 

-0.10 


32.8948 

3573.927 

16.66 


0.11 

0.272 

0.0039 

0.0038 


mm 

1.02 

7992 

0.16 


49.8658 

3624.125 

16.75 



0.477 

0.0052 

0.0005 


mm 

1.06 

8612 

0.02 


27.0844 

3607.949 

16.27 



0.462 

0.0079 

0.0054 


mm 

2.82 

12973 

-0.07 


30.8419 

3590.014 

16.69 


0.15 

0.511 

0.0045 

0.0060 


mm 

1.63 

12987 

0.02 


25.9121 

3584.179 

17.41 


0.13 

0.363 

0.0034 

0.0065 


2 

1.08 

14083 

0.01 


29.9824 

3583.922 

17.51 


gilETil 

0.657 

0.0054 

0.0045 


2 

1.03 

17017 

-0.10 


25.7548 

3637.727 

16.29 


0.11 

0.519 

0.0041 

0.0007 


mm 

1.10 

18037 

0.00 


31.5140 

3563.204 

17.00 


gflETil 

0.573 

0.0047 

0.0063 


mm 

1.46 

19353 

0.02 


22.2690 

3587.215 

17.14 


0.34 

0.448 

0.0077 

0.0058 


mm 

1.24 

20313 

0.04 


27.7087 

3602.991 

16.40 


0.59 

0.496 

0.0012 

0.0031 


mm 

1.38 

23350 

0.07 


28.3841 

3623.853 

17.40 


0.63 

0.369 

0.0029 

0.0011 


2 

1.10 


Category 7: list of 2 peculiar EBs that exhibit variations between eclipses. 


] Catalog Properties 

Analytic Model Parameters 

Fit Statistics j 

ID 

(V-I) 


p 

to 

(I) 

Ah 

Ah 

^>2 

©1 

02 

fa.I 

iEM 

XcJv 

343 

0.13 


33.5508 

3572.393 

17.35 


0.17 

0.694 

0.0124 

0.0087 

BlETtl 

1 

1.03 

4458 

-0.07 


39.5436 

3581.850 

16.14 


[jgg] 

0.501 

0.0029 

0.0007 


2 

1.41 
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Table 1 (cont.): 

Category 8: list of 130 detached EBs with well-defined eclipse parameters. 


Catalog Properties 


Analytic 

Model Parameters 



Fit Statistics 

ID 

(V-I) 

Mi 

P 

to 

(I) 

Ah 

Ah 

$2 

01 

©2 

fa ,/ 



91 

0.06 

411 

24.8099 

3619.170 

16.80 

gjyjii 


0.459 

0.0080 

0.0114 



1.02 

170 

-0.03 

426 

26.3720 

3566.276 

17.10 



0.155 

0.0036 

0.0066 



1.05 

784 

0.00 

444 

44.1196 

3633.801 

16.72 



0.115 

0.0026 

0.0035 


mm 

1.02 

866 

0.06 

424 

28.0132 

3598.047* 

16.47 



0.496 

0.0046 

0.0101 



1.04 

1056 

0.01 

444 

30.6713 

3647.862* 

17.40 



0.503 

0.0037 

0.0065 



1.04 

1530 

0.06 

465 

42.1989 

3632.166 

17.42 



0.655 

0.0044 

0.0033 


MM 

1.02 

1968 

0.13 

911 

44.8940 

3627.360 

17.33 



0.437 

0.0035 

0.0038 



1.01 

2142 

0.01 

457 

27.7692 

3595.213 

16.57 



0.162 

0.0047 

0.0050 



1.05 

2277 

0.01 

465 

36.5936 

3616.509 

17.08 



0.473 

0.0039 

0.0041 


BB 

1.00 

2708 

0.04 

465 

44.5975 

3633.684 

17.48 



0.422 

0.0066 

0.0031 


1 

1.03 

2780 

-0.12 

446 

27.0438 

3628.748* 

17.03 



0.375 

0.0024 

0.0088 



1.02 

3082 

-0.06 

876 

45.5327 

3589.237 

16.71 



0.710 

0.0024 

0.0024 


mm 

1.00 

3177 

-0.12 

434 

20.2216 

3585.258 

16.70 



0.455 

0.0045 

0.0106 


1 

1.04 

3388 

-0.08 

448 

44.9697 

3588.586* 

16.53 



0.168 

0.0019 

0.0049 


1 

1.00 

3557 

0.00 

412 

23.8734 

3592.035 

17.03 



0.478 

0.0048 

0.0062 


1 

1.01 

4031 

-0.10 

448 

32.5105 

3577.358 

16.46 



0.821 

0.0045 

0.0035 



1.04 

4399 

0.03 

447 

22.9067 

3624.060* 

16.76 



0.644 

0.0068 

0.0074 


mm 

1.03 

4419 

0.07 

434 

42.0950 

3575.980 

16.09 



0.415 

0.0032 

0.0046 


1 

1.02 

4721 

-0.05 

456 

21.8793 

3592.692 

17.05 



0.511 

0.0059 

0.0073 


mm 

1.06 

4737 

0.07 

440 

31.5717 

3612.900* 

17.51 



0.644 

0.0036 

0.0080 


1 

1.00 

4804 

-0.02 

429 

23.0375 

3601.271 

16.23 



0.463 

0.0046 

0.0048 


1 

1.00 

4837 

-0.02 

445 

26.8055 

3608.319* 

17.56 



0.576 

0.0050 

0.0043 


mm 

1.01 

5145 

0.00 

445 

47.6341 

3684.851 

17.55 



0.173 

0.0036 

0.0033 



1.00 

5153 

-0.04 

456 

24.4530 

3582.856* 

16.28 



0.829 

0.0039 

0.0031 


1 

1.02 

5195 

-0.06 

445 

39.4245 

3564.482 

17.19 



0.439 

0.0064 

0.0047 


mm 

1.06 

5492 

-0.04 

429 

32.9664 

3602.498 

16.27 



0.280 

0.0042 

0.0085 


1 

1.03 

5965 

-0.05 

439 

20.3374 

3616.353 

17.54 



0.684 

0.0041 

0.0059 


mm 

1.02 

6187 

0.00 

477 

29.8706 

3612.492* 

16.91 



0.537 

0.0042 

0.0085 


1 

1.08 

6555 

-0.02 

477 

29.0548 

3591.351 

17.19 



0.657 

0.0071 

0.0038 


1 

1.00 

6996 

-0.04 

476 

29.9285 

3596.250* 

16.54 



0.445 

0.0041 

0.0096 


mm 

1.05 

7380 

0.02 

468 

31.4810 

3597.756* 

17.23 



0.564 

0.0036 

0.0045 


mm 

1.00 

7560 

-0.01 

418 

26.7461 

3611.788* 

16.49 



0.575 

0.0036 

0.0032 


mm 

1.09 

7565 

0.13 

477 

29.0955 

3601.070 

17.45 



0.667 

0.0103 

0.0051 



1.01 

7935 

-0.01 

404 

24.0615 

3609.910 

17.39 



0.449 

0.0045 

0.0051 


1 

1.05 

7975 

-0.05 

477 

20.1226 

3615.517* 

16.86 



0.607 

0.0031 

0.0072 


mm 

1.03 

8543 

0.00 

451 

25.3380 

3621.313 

17.07 



0.414 

0.0067 

0.0057 


1 

1.05 

8559 

-0.06 

468 

42.2643 

3631.168* 

17.54 



0.239 

0.0019 

0.0045 


mm 

1.04 

8783 

0.08 

458 

48.7243 

3622.458 

17.44 



0.170 

0.0041 

0.0041 


1 

1.06 

8903 

-0.05 

463 

25.6605 

3637.923* 

16.49 



0.593 

0.0049 

0.0050 


1 

0.99 

8993 

-0.09 

891 

29.4599 

3599.445* 

17.11 



0.617 

0.0024 

0.0073 


mm 

0.99 

9159 

0.05 

466 

48.1575 

3695.306 

17.34 



0.622 

0.0038 

0.0054 



1.00 

9386 

-0.07 

559 

29.8063 

3619.572* 

16.61 



0.682 

0.0030 

0.0062 


1 

1.00 

9429 

0.11 

476 

20.6681 

3605.574* 

17.15 



0.532 

0.0075 

0.0126 


1 

1.08 

9441 

0.05 

476 

20.8997 

3596.643 

17.45 



0.432 

0.0058 

0.0075 



1.00 

9953 

-0.10 

490 

45.5783 

3729.995* 

17.05 



0.140 

0.0025 

0.0032 


mm 

1.05 

10096 

0.05 

431 

35.2251 

3628.920 

17.58 



0.317 

0.0049 

0.0060 


1 

1.01 

10422 

-0.11 

612 

21.1694 

3586.720* 

16.91 



0.841 

0.0032 

0.0045 


1 

1.02 

10575 

-0.05 

559 

29.1436 

3592.083 

17.14 



0.546 

0.0057 

0.0049 


mm 

1.04 

10953 

-0.12 

616 

42.2762 

3610.382 

16.31 



0.265 

0.0059 

0.0044 


1 

1.08 

11252 

-0.04 

490 

22.2402 

3592.677 

17.19 



0.496 

0.0067 

0.0052 


1 

1.02 

11299 

0.07 

1049 

39.5635 

3584.911* 

17.22 



0.214 

0.0055 

0.0063 


1 

1.01 

11526 

-0.01 

477 

49.2837 

3595.345 

16.67 



0.418 

0.0035 

0.0041 



1.09 

11538 

-0.17 

477 

33.7813 

3600.922* 

16.28 



0.504 

0.0017 

0.0092 


mm 

1.02 

11636 

0.08 

474 

31.1510 

3579.485 

17.59 



0.612 

0.0054 

0.0069 


1 

1.09 

11907 

-0.04 

440 

45.4195 

3629.747 

17.20 



0.589 

0.0031 

0.0047 



1.00 

12170 

-0.04 

452 

29.1744 

3586.509 

17.60 



0.439 

0.0060 

0.0045 



1.01 

12179 

-0.01 

477 

25.1020 

3555.863* 

16.84 



0.485 

0.0032 

0.0098 



0.99 

12384 

0.12 

476 

25.1786 

3591.906 

17.12 



0.251 

0.0116 

0.0096 


mm 

1.08 

12454 

-0.14 

476 

30.1523 

3591.689 

17.01 



0.649 

0.0043 

0.0029 


E3I 

1.04 

12832 

-0.08 

493 

20.1668 

3551.653* 

16.37 



0.425 

0.0047 

0.0049 


— 

1.04 

13177 

-0.01 

477 

40.5158 

3542.700 

17.11 



0.813 

0.0034 

0.0058 


mm 

1.00 

13260 

0.00 

450 

24.2957 

3593.162 

17.21 



0.447 

0.0095 

0.0066 



1.16 

13390 

0.10 

477 

35.7858 

3571.607 

17.59 



0.630 

0.0026 

0.0057 


1 

1.01 

13418 

-0.11 

465 

26.9111 

3587.289* 

16.94 



0.907 

0.0025 

0.0040 



1.03 

13441 

-0.04 

493 

30.6531 

3593.354 

17.26 



0.618 

0.0042 

0.0044 


1 

1.03 


*Epoch of primary eclipse minimum t 0 appropriately adjusted to ensure A/i > A 1 2 ■ 
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Table 1 (cont.): 

Category 8 (cont.): list of 130 detached EBs with well-defined eclipse parameters. 


Catalog Properties 


Analytic 

Model Parameters 



Fit Statistics 

ID 

(V-I) 


P 

to 

(I) 

A/i 

A I 2 

<|>2 

©i 

©2 

f<7,I 

Hi i 


13482 

-0.08 


21.5789 

3548.902* 

17.32 

amim 

mgKi 

0.659 

0.0050 

0.0058 

1.10 

mm 

0.99 

13491 

0.10 


21.3014 

3578.036* 

16.82 



0.510 

0.0069 

0.0069 

mi 

2 

1.07 

13726 

0.10 


20.9344 

3554.051 

17.37 



0.679 

0.0050 

0.0069 


2 

1.03 

13867 

-0.02 


21.7457 

3604.874 

17.33 



0.493 

0.0047 

0.0051 


1 

1.00 

14171 

-0.02 


21.7291 

3621.792 

16.95 



0.528 

0.0043 

0.0051 


MM 

1.02 

14360 

0.00 


30.2612 

3651.363* 

16.20 



0.595 

0.0038 

0.0072 


1 

0.99 

14895 

0.06 


34.4641 

3574.345* 

17.33 



0.186 

0.0040 

0.0034 



1.01 

15235 

-0.10 


28.9500 

3553.091* 

17.17 



0.317 

0.0045 

0.0047 


MM 

1.05 

15244 

0.07 


21.8435 

3622.895 

17.41 



0.299 

0.0064 

0.0071 


2 

1.02 

15380 

-0.09 


27.6171 

3624.679* 

16.79 



0.802 

0.0028 

0.0045 



1.05 

15788 

-0.06 


29.0140 

3561.929* 

17.04 



0.246 

0.0031 

0.0033 

1.27 

MM 

1.04 

15979 

-0.12 


28.1428 

3635.572 

16.27 



0.637 

0.0077 

0.0032 

area 

1 

1.08 

16026 

0.00 


30.8480 

3566.099 

16.40 



0.661 

0.0038 

0.0103 



0.87 

16126 

0.02 


22.2111 

3609.525 

17.49 



0.234 

0.0038 

0.0075 


1 

0.99 

16350 

-0.10 


30.9801 

3603.346* 

17.09 



0.746 

0.0027 

0.0046 



1.02 

16399 

0.01 


39.1458 

3653.057* 

17.45 



0.849 

0.0020 

0.0039 



1.00 

16418 

0.03 


20.4164 

3586.945 

17.08 



0.619 

0.0070 

0.0079 

1.05 


1.03 

16711 

-0.13 


26.5576 

3632.458* 

16.25 



0.833 

0.0027 

0.0056 

BiiEl 

1 

1.04 

16964 

-0.10 


23.5422 

3603.149 

16.44 



0.846 

0.0031 

0.0052 



1.03 

17067 

-0.14 


28.1151 

3592.517 

16.57 



0.666 

0.0033 

0.0070 



1.02 

17316 

-0.14 


26.7747 

3611.392 

16.68 



0.775 

0.0056 

0.0027 


mm 

1.00 

17361 

0.03 


20.9318 

3562.002 

17.48 



0.551 

0.0111 

0.0065 



0.98 

17539 

0.06 


23.9094 

3562.789* 

17.51 



0.471 

0.0053 

0.0082 


2 

0.99 

17569 

-0.08 


22.2019 

3591.019 

16.91 



0.523 

0.0049 

0.0035 


MM 

1.12 

17750 

-0.11 


20.7086 

3586.855 

16.98 



0.750 

0.0066 

0.0072 


1 

1.01 

17784 

-0.09 


25.2613 

3575.400 

17.11 



0.280 

0.0038 

0.0031 


1 

1.01 

17822 

0.07 


44.2706 

3641.930 

16.75 



0.579 

0.0072 

0.0048 

1.09 

1 

1.01 

18237 

-0.10 


20.3285 

3616.375 

16.03 



0.481 

0.0105 

0.0092 

BIEE8 

1 

1.14 

18582 

0.02 


41.9856 

3608.318* 

16.60 



0.539 

0.0047 

0.0065 


MM 

1.06 

18659 

0.04 


21.4449 

3594.435 

17.14 



0.665 

0.0079 

0.0048 


1 

1.03 

18813 

0.11 


39.6989 

3627.894* 

16.92 



0.498 

0.0041 

0.0078 



1.02 

18824 

-0.08 


33.5726 

3576.478 

16.54 



0.439 

0.0044 

0.0064 


1 

1.03 

18839 

-0.04 


48.2133 

3619.336* 

17.28 



0.269 

0.0037 

0.0036 


1 

0.91 

18859 

0.00 


24.5592 

3624.022 

17.01 



0.505 

0.0038 

0.0096 


0 

1.00 

18869 

-0.08 


42.5786 

3628.587 

16.90 



0.648 

0.0039 

0.0023 


0 

1.00 

19083 

0.02 


22.2719 

3600.521 

16.85 



0.517 

0.0062 

0.0062 


MM 

1.09 

19230 

0.04 


37.4005 

3569.116 

16.89 



0.280 

0.0040 

0.0041 


0 

0.98 

19792 

-0.10 


29.6438 

3640.465 

16.00 



0.519 

0.0047 

0.0042 

1.52 

0 

1.01 

19840 

0.02 


30.9709 

3617.985 

16.15 



0.449 

0.0080 

0.0046 

BlfcfrM 

1 

1.04 

20309 

-0.02 


48.2396 

3694.899 

17.60 



0.443 

0.0034 

0.0039 


0 

1.01 

20459 

-0.14 


30.7408 

3642.410* 

17.22 



0.590 

0.0033 

0.0039 


MM 

1.02 

20522 

0.07 


27.2575 

3619.621 

17.40 



0.457 

0.0058 

0.0073 


1 

1.06 

20590 

0.06 


40.0925 

3618.364* 

17.32 



0.563 

0.0063 

0.0058 

1.07 

0 

1.11 

20646 

0.04 


26.0462 

3620.172 

16.91 



0.656 

0.0068 

0.0066 

Bills* 

0 

1.02 

20746 

0.08 


28.3719 

3614.200* 

16.69 



0.575 

0.0045 

0.0090 


0 

1.02 

21059 

0.07 


25.7904 

3611.136* 

17.56 



0.674 

0.0043 

0.0046 


0 

0.98 

21518 

0.02 


21.4601 

3629.580 

17.10 



0.443 

0.0074 

0.0077 


MM 

1.08 

21621 

0.04 


38.5314 

3576.469* 

17.20 



0.341 

0.0032 

0.0075 


2 

0.96 

21881 

0.04 


21.0056 

3588.139* 

16.80 



0.354 

0.0050 

0.0073 

1.10 

0 

1.05 

22082 

-0.06 


33.1149 

3586.199 

16.76 



0.341 

0.0056 

0.0027 


0 

1.04 

22553 

0.01 


22.8430 

3598.817 

17.31 



0.804 

0.0042 

0.0061 


2 

0.96 

22691 

0.00 


31.8865 

3579.266* 

16.89 



0.235 

0.0018 

0.0068 


0 

0.98 

22713 

-0.09 


33.3752 

3617.549 

17.00 



0.679 

0.0037 

0.0035 


1 

1.02 

22764 

0.00 


29.7044 

3572.606* 

16.63 



0.155 

0.0026 

0.0050 

1.13 

0 

1.03 

23088 

0.09 


22.4312 

3598.919 

17.43 



0.495 

0.0065 

0.0070 

■EH 

0 

1.07 

23101 

0.08 


46.3819 

3616.728* 

16.79 



0.196 

0.0014 

0.0094 


1 

1.06 

23368 

0.03 


22.9011 

3629.347* 

17.29 



0.646 

0.0037 

0.0056 


1 

0.99 

23773 

0.06 


35.2308 

3570.982* 

16.75 



0.372 

0.0040 

0.0055 


0 

1.08 

24195 

0.13 


36.6221 

3639.504 

17.33 



0.235 

0.0041 

0.0047 

1.05 

0 

1.02 

24580 

-0.06 


31.1688 

3608.707 

16.50 



0.456 

0.0043 

0.0027 

QQ 

1 

1.01 

24818 

-0.04 


21.5495 

3586.096 

17.54 



0.778 

0.0050 

0.0027 


0 

1.02 

24858 

0.13 


45.0560 

3623.442 

17.14 



0.937 

0.0029 

0.0033 


1 

1.03 

25297 

0.04 


34.3872 

3647.584 

17.44 



0.225 

0.0053 

0.0055 


0 

1.03 

25578 

-0.11 


21.7054 

3582.818 

16.08 



0.298 

0.0060 

0.0097 

1.34 

MM 

1.01 

26109 

-0.04 


34.5497 

3622.492 

16.55 



0.864 

0.0035 

0.0041 

■IM 

1 

0.99 


*Epoch of primary eclipse minimum t 0 appropriately adjusted to ensure A/i > A 1 2 ■ 
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Table 2: Physical model properties and statistics for the 130 detached EBs in the well-defined sample. We list the 
OGLE-III LMC EB identification number and the nine physical model properties: orbital period P (days), epoch of primary 
eclipse t 0 (JD — 2450000), primary and secondary component masses M p and M a (Mg), age r (Myr), inclination i (°), 
eccentricity e, argument of periastron uj (°), and dust extinction Ai (mag). We then list other physical properties including the 
mass ratio q = M 2 /M 1 = min{M p , M s }/max{M p ,M s }, orbital separation a (R©), stellar radii R p and R s (R©), and effective 
temperatures T p and T s (K). Finally, we list the fit statistics including the photometric correction factors and f a ,v, number 
of data points Mi and Mv, and number of data points we clipped M c ,i and M c ,v in the I-band and V-band, respectively. 


Independent Physical Model Properties Dependent Physical Properties | Fit Statistics 


12,500 11,900 


ID 

P 


91 

24.8098 

3619.171 


26.3719 

3566.280 

784 

44.1198 

3633.792 


28.0132 

3598.058 


30.6714 

3647.865 


42.1990 

3632.165 

1968 

44.8939 

3627.357 


27.7691 

3595.207 


36.5937 

3616.510 


44.5972 

3633.673 


27.0440 

3628.753 


45.5329 

3589.236 


20.2215 

3585.256 


44.9699 

3588.605 


23.8734 

3592.036 


32.5105 

3577.365 


22.9066 

3624.073 

4419 

42.0952 

3575.980 


21.8791 

3592.698 

4737 

31.5717 

3612.891 

gain 

23.0376 

3601.269 

4837 

26.8054 

3608.320 

5145 

47.6339 

3684.852 


24.4530 

3582.859 

5195 

39.4239 

3564.495 

■am 

32.9666 

3602.509 

5965 

20.3374 

3616.352 

laiai 

29.8703 

3612.495 

6555 

29.0549 

3591.353 


29.9282 

3596.246 


31.4810 

3597.752 


26.7462 

3611.776 

7565 

29.0954 

3601.066 


24.0616 

3609.906 


20.1227 

3615.519 

8543 

25.3381 

3621.308 


42.2644 

3631.167 


48.7240 

3622.451 


25.6604 

3637.910 

8993 

29.4599 

3599.449 


48.1568 

3695.292 


29.8063 

3619.568 


20.6681 

3605.551 

9441 

20.8997 

3596.647 


45.5780 

3730.003 


35.2251 

3628.935 


21.1695 

3586.712 


29.1431 

3592.096 


42.2762 

3610.406 


22.2402 

3592.679 

11299 

39.5635 

3584.917 

Hfe'ii 

49.2822 

3595.336 

11538 

33.7813 

3600.922 


31.1511 

3579.505 


45.4194 

3629.746 


29.1745 

3586.504 

12179 

25.1022 

3555.832 


25.1787 

3591.907 


30.1522 

3591.701 

12832 

20.1669 

3551.657 

13177 

40.5163 

3542.688 

■KWiill 

24.2959 

3593.154 


13,400 13,300 


15,400 12,700 


15,100 11,800 


13,200 12,500 


17,200 12,300 


11,100 


19,900 12,300 


20,500 13,700 


14,500 15,700 


22,600 23,100 































































































































































































































































































































































































